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conveying  any  rights  or  permission  to  manufacture  use,  or  sell  any  patented  invention 
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1.  INTRODUCTION 


This  report  describes  the  structure  and  use  of  the  program  PARCS  (Plume 
Attenuated  Radar  Cross  Section).  This  is  a CDC  FORTRAN  IV  code  that  calculates 
the  coherent,  incoherent,  and  overdense  surface  radar  cross  sections  of  a rocket 
plume.  The  calculations  include  attenuation,  local  index  of  refraction,  doppler 
shift,  and  range  cell.  The  program  accepts  plume  data  directly  from  the  AeroChem 
LAPP  code  but  may  be  interfaced  with  other  sources  of  plume  definition. 

The  program  was  developed  for  the  Air  Force  Rocket  Propulsion  Laboratory 
under  contract  No.  F04611  -75-C-0021,  monitored  by  Capt.  W.  Rothschild  and 
Lt.  R.  Sperlein.  The  Technical  Monitor  was  Dr.  R.  Fante  (AFGL). 

The  RCS  model  was  developed  at  Aerocfyne  by  Drs.  D.  Mann  and  J.  Draper/^*  ^ 
This  report  summarizes  the  model  equations  and  gives  an  outline  of  the  program. 


1.  R.  R.  Mikatarian,  C.  J.  Kau  and  II.  S.  Pergament,  "A  Fast  Computer  Program 
for  Nonequilibrium  Rocket  Plume  Predictions,  "AeroChem  Research  Laboratories, 
Inc.,  August  1972,  AFRPL-TR-72-94. 

2.  1).  Mann,  J.  Draper,  and  J.  Rickman,  '’Rocket  Plume  Radar  Cross  Sections:  Theory 
and  Data."  JANNAF  9th  PLUME  TECHNOLOGY  MEETING. 

3.  JANNAF  Rocket  Exhaust  Plume  Technology  Hamfcook,"  Plume  RCS,"  J.S.  Draper 
and  1).  Mann,  soon  to  be  published. 
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2.  DESCRIPTION  OF  THE  MODEL 


2 . 1 Table  of  Symbols 


Text 

Program 

a 

ATT 

Attenuation  coefficient  (m~S 

c 

C 

Velocity  of  light  (2.998  x 10®  m/sec) 

f 

RFREQ 

Radar  frequency  (Hz) 

4f 

Doppler  frequency  shift  (Hz) 

4f 

D 

DOPFRQ 

Radar  frequency  resolution  (Hz) 

K 

Radar  wave  vector  (m~*) 

L 

FL 

Maximum  plume  length  intersected  by  radar  beam 

A 

N 

(RJET) 

Normal  to  overdense  surface 

N 

Complex  Index  of  refraction 

N 

XI ND 

Real  part  of  index  of  refraction 

ft 

YIND 

Negative  of  imaginary  part  of  Index  of  refraction 

Nei,k 

PLE  (I,  K) 

Plume  electron  density  (cm  ) 

Ne 

ZNE 

_3 

Interpolated  electron  density  (cm  ) 

P 

PHASE 

Coherent  cross  section  phase  (radians) 

fj 

RJET 

Jet  radius  at  matched  pressure  (m) 

rU 

PLR  <1,  J) 

Radial  location  of  plume  data  point  (m/jet  radius) 

rodj 

PLROVD(l.J) 

Radius  of  OVERDENSE  surface  at  z^  (m/jet  radius) 

R 

RANGE 

Radar  source  to  plume  range  (m) 

s 

S 

Ray  S coordinate 

S . S 
min.  max 

SMIN,  SMAX 

Ray  integration  limits 

V 

PLVC0.J) 

Plume  gas  axial  velocity  (m/sec) 

VD 

DOPVEL 

Radar  velocity  resolution 

W 

B WIDTH 

Radar  beam  width  (m) 

xozo 

XO.ZO 

Ray  Intersection  with  x.  z - plane 

xLlm 

XUM 

x coordinate  Integration  limit 

z00  * Urn 

ZOO,  ZUM 

z coordinate  Integration  limit 
2 

I 

\ 

\ 


Text 

Program 

Forward  scattering  angle  (radians) 

y 

GAM 

Mean  fluctuation  factor 

ds 

STEP 

Step  size  along  s-axis  (jet  radii) 

O 

X 

DXO 

Step  size  along  x-axis  (jet  radii) 

4zo 

DZO 

Step  size  along  z-axis  (jet  radii) 

e 

THETA D 

Radar  aspect  angle  (DEG)  10°  nose  on  viewing) 

«>w 

FMUOD(l), 

FMU0D(2) 

_2 

Complex  overdense  surface  return  amplitude  (m  ) 

"ij 

PLFC(I.J) 

Plume  collision  frequency 

Tt 

PI 

3. 141592654 

P 

Overdense  surface  reflection  coefficient  amplitude 

a 

SIG 

2 

Total  ranar  cross  section  (RCS)  (m  ) 

^coh 

SCOH 

Coherent  RCS  (m2) 

^inc 

SING 

incoherent  RCS  (m“) 

rod 

SOD 

Overdense  RCS  (m'l 

TTI LDA  ( 1 > , 
TT1LI)A(2) 

_•> 

Coherent  RCS  complex  amplitude  (m  ) 

* 

Angle  between  overdense  surlnce  x,z  plane 
projection  and  the  z-axis 

CJ 

OMKG 

Radar  frequency  (radians) 
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2.2  Coordinate  Systems 

The  rocket  plume  coordinate  system  is  shown  in  Figure  1.  The  plume  is 
assumed  to  be  symmetric  about  the  a -axis.  The  x,  y plane  is  parallel  to  the 
nozzle  exit  plane,  bu*.  displaced  by  an  amount  dependent  on  the  interval  over 
which  the  plume  data  is  given.  The  coordinate  values  are  in  units  of  jet  radius 
at  matched  pressure. 

The  plume  data  is  given  at  discrete  points  r^  by  the  AeroChem  LAPP 

code!1)  The  plume  properties  given  are  electron  density  N , electron-neutral 

0 

collision  frequency  v and  axial  gas  velocity  v. 

Volume  and  surface  integrations  are  performed  in  the  skew  coordinate 
system  obtained  by  rotating  the  plume  system  y-axis  so  that  it  lies  parallel 
to  the  radar  line-of-sight.  The  rotated  axis  is  called  the  s-axis,  and  it  can  be 
seen  in  Figure  1 that  y = -s  sin  0 , where  0 is  the  radar  aspect  angle. 

Volume  integrations  are  performed  over  a matrix  of  rays  parallel  to  the 
s-axis,  and  a volume  dv  = dx  dz  = cbt  dz  (-ds  sin  0 ) is  associated  with  each 
integration  step. 

2.3  Incoherent  RCS 

The  underdense  Incoherent  cross  section  is  calculated  from  the  formulas: 


r 


\ 


i 


¥ 


ds  • 


| 
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Figure  1 - Coordinate  System 

llie  integration  limits  are  given  by: 
x.  . PL  KM  AX 

Llm 

?-00  - -PLKMAX/tan  6 

? j im  P LRM AX/tan  9 + I.  , 


where  PLRMAX  is  the  maximum  plume  radius,  L is  the  length  of  the  plume,  and 
9 is  the  aspect  angle. 


Other  definitions  required  are: 


The  quantity  "9  is  the  (nondimensional)  radar  cross  section  per  unit  area 

along  the  ray  that  passes  through  (xQ,  0,  zQ).  The  limits  of  integration  are 

obtained  by  Intersecting  the  ray  with  a plume  defined  as  a set  of  cylindrical 

slabs  (see  Figure  2).  The  integrand  In  the  9 calculation  is  the  product  of  a 

nondimensional  attenuation  factor  and  the  cross  section  per  unit  volume 
-1  -1  3 

do  / dfl  (m  sr  ).  The  function  #(m  ) is  the  Kolmogorov  correlation 
function,  and  a (m  *)  is  the  attenuation  coefficient  per  unit  length. 

Other  definitions  required  are: 

•15 

rfi  = electron  radius  = 2.3178x10  (m) 

y - mean  fluctuation  factor 

turbulence  scale  - r^  (m)  (Jet  radius  at  matched  pressure) 

c - velocity  of  light  * 2.998  x 108  (m/sec) 

u - signal  frequency  *-  2rrf  (sec  l) 

2 

N = square  of  index  of  refraction 

Nc  » critical  density  = 3. 14  x 10’10  (cm  *) 

0^  m forward  scatter  1/2  angle  * n 

The  mean  fluctuation  factor,  y.  Is  the  ratio  of  the  fluctuating  to  the  mean 
electron  density,  lyplcal  values  are  between  0.  5 and  1.  0 
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2.4  Coherent  RCS 

The  underdense  coherent  cross  section  is  calculated  from  the  formulas: 


ZLim  XLim 


where 


47r  / J 


0 0 


yl1  + (tj-)2 


— exp 


-2  ^ arjds'J  (cos  P-i  sin  P)  • r^  sin#  ds 
Smin  / (ID 


a 

p ’ / rjdB'  + 2 f r)(smin  + z0  005  9 ) • 


This  calculation  is  complex  to  include  both  amplitude  and  phase.  The  quantity  |f| 
has  the  units  of  m *.  All  other  symbols  have  been  defined  in  Subsection  2.2.  The 
quantity  (s  + zQ  cos  g ) is  the  distance  from  the  edge  of  the  plume  (smjn  in 
Figure  2)  to  a fixed  plane  passing  through  the  origin  normal  to  the  integration  ray. 
Those  rays  that  pass  through  the  plume  rocket  interface  (Figure  2)  are  made  to  give 
a zero  coherent  and  overdense  return  because  the  plume  face  is  obscured  by  the 


2.5  Overdense  RCS 

The  integrations  in  Subsections  2.3  and  2.4  with  respect  to  s are  terminated  when 
an  overdense  surface  is  reached.  This  is  defined  to  be  the  locus  of  points  where. 


N,  =Nod  -Nc  [>  '{irf]  [l  *(£)]" 


sin20 
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The  radar  cross  section  contribution  from  reflection  by  the  overdense 
surface  is  given  by: 


od  ' 7T 


ZLim  XLim 


/ / 2"(xo’zo)rj2  8indxi 


odzo 


(14) 


with 


/V  As 

'v  • K 


P = ^ ^1/2  exp  f -2  J ar^ds  J (cos  P - i sin  P) 


min 


(15) 


1 - N 
~ 1 + N 


(N  e complex  index  of  refraction) 


(16) 


fi  = 0.  8 r^  sin  6 


radius  of  overdense  surface 

(17) 


A A 

where  K is  the  unit  wave  vector,  N is  the  normal  to  the  overdense  surface,  p is  the 
reflection  coefficient  amplitude,  and  £ is  a surface  roughness  correction.  The 
surface  roughness  becomes  important  when  the  mean  roughness  (h)  is  on  fie 
order  of  a radar  wavelength.  Ir  the  above  formulation,  h - 0.  4 and  surface 
roughness  should  be  included  for  kh  > 0.  5. 

The  radii  of  the  overdense  surface  r . at  each  axial  station  i,  are 

od  1 

determined  from  the  overdense  condition  1 

i Nc  [>  ♦ (-C-)2]  ( 1 * 2-fc)'1  “l"2s  • 

(18) 
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The  intermediate  values  rQ(j  (z)  are  linearly  interpolated  from  the  plume  axial 
station  values  rQ(j  and  rQ(j  , where, 

Zj  _ j C Z < Zj  • 

The  normal  to  the  overdense  surface  in  the  x z-plane  is 

N(z)  = cos  , 0,  -sin  \p  , (19) 

where 

r°d{+ 1 rodj 

tan  = , with  z.  < z < z.  * 1 . (20) 

zi+l  Zi  11 

Hie  three  dimensional  normal  at  an  arbitrary  point  (x,  y,  z)  is  obtained  by 
rotation  about  the  z-axls. 


N(x,y,  z)  - (cos  \l/  cob  0,  cos  ^sin  0,  - sin  \Jf)  , 


for 


cos  0 


4 


2 2 

x + y 


sin  0 


I 2 2 

+ y 


The  radar  wave  vector  dot  product  with  the  surface  normal  is  given  by: 


N • K -sin  9 cos  \l>  sin  0 - cos  9 sin  , 


(21) 

(22) 


(23) 


where 


K (0,  -sin  9,  cos  9) 


(24) 


I 
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2.  6 Range  Cell 

The  radar  resolution  is  limited  by  the  range  cell  length  (/)  and  beam  width 
(w)  as  shown  in  Figure  3.  For  small  aspect  angles  (0)  the  resolution  is  limited 
by  the  range  cell  length,  and  the  maximum  plume  length  intersected  by  the  radar 
beam  is 


L = 


w 

sin  $ 


0 < 


6 < tan 


(25, 


Figure  3 - Finite  Radar  Beam  Incident  on  Cylindrical  Plume 
(for  w = / -=  104  6 3 tan-1  (w // )) 


11 


For  large  aspect  angles  the  maximum  plume  length  intersected  is 
determined  by  the  range  cell  length. 

L ■ co/t  ■ ta"'1(r)  s « i f ■ <26> 

The  large  and  small  aspect  angle  equations  are  equal  for  B = tan-1(a//)  as 
shown  in  Figure  3.  The  beam  width  is  determined  from  the  range  (R)  and 
beam  divergence  (6)  according  to 

w = R 6 . (27) 

The  cross  sections  are  summed  and  printed  for  each  range  cell  when  the 
integration  intervals  in  the  axial  direction  (z  in  Figure  3)  are  multiples  of  the 
plume  intersection  length  L.  The  low  aspect  angle  case  (0  < 0Q)  Is  not  range 
cell  resolved,  but  appears  in  a single  range  cell  (see  Subsection  2.8). 

2.  7 Doppler  Shift 

The  incoherent  cross  section  is  determined  as  a function  of  doppler 
frequency  shift.  The  shift  is  due  to  the  local  gas  velocity  (\)  relative  to  the 
missile,  and  is  given  by: 

df(Hz)  = 2ttK  • V , (28) 

where  K is  the  radar  wave  vector.  The  contributions  to  the  incoherent  cross 
section  at  each  point  along  the  integration  ray  are  put  into  a Doppler  bin  deter- 
mined by  the  local  gas  velocity. 

The  velocity  resolution  of  the  radar. 

4fD  e Radar  frequency  resolution  (29) 


D 


4_Id 

2*rK 
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and  the  maximum  gas  velocity  are  used  to  define  the  velocity  width  of  the 
Doppler  bins  (Vn)  and  the  number  of  bins.  As  the  integration  along  the  ray 
proceeds  the  local  gas  velocity  V is  projected  along  the  radar  line-of-s'ght  K 
and  a Doppler  bin  N is  selected  such  that 

(N-l)  VD  < V • K < N VD  . (30) 

The  local  incoherent  cross  section  contribution  is  stored  in  Doppler  bin  N 
and  is  taken  to  be  the  integral  along  the  ray  from  the  preceding  integration 
point. 

2.  8 Low  Aspect  Angle 

The  cross  sections  for  aspect  angles  in  the  range 

j$j  < $ (0Q  = low  aspect  angle  definition  from  Input  data)  , 

(31) 

are  determined  by  integrating  over  a matrix  of  rays  defined  in  the  x,  y-plane. 

(The  rays  are  defined  in  the  x,  z-plane  for  higher  aspect  angles. ) The  coherent 
and  incoherent  cross  sections  are  determined  from  the  equations  of  Subsections 
2.  3 and  2.  4 with  the  substitutions 

and 

sin  $ ds  — ►(  t"  ■ z)  ds  , (32) 

where  T is  a unit  vector  tangent  to  the  ray  and  s is  the  path  length. 
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The  Integration  ray  bends  as  it 


steps  through  the  plume  according  to 


where 


t(P') 


(33) 


£-(***)  8-*  i(fv|Bi)B[ 

♦*  i(£v  £».)». 


and 


U = X t = 

\ VN  ! 


1 

PN 


(34) 


- (t  gN\  1 ^ 

or  [ rdr  lz  fr  ) lr  J r 

. (t  5N;+  t 0N\t  1 £_  (t  dN  aN\  * 

\ror  ‘zazJ^J2  ‘zarj1^  ^ 


(35) 


where  r is  a unit  vector  along  the  cylindrical  radius  from  the  plume  axis, 

<p  is  perpendicular  to  r such  that  $ x r * z and  PN  is  the  real  part  of  the 
gradient  of  the  index  of  refraction. 


The  gradient  of  the  index  of  refraction  is  determined  from: 


14 


(39) 


The  initial  value  of  t is  determined  from 

dR  = tr‘r  + t^£+  't^z  ds  = (-sin  $ y + cos  9 z)  ds 


to  be 

(40) 

(41) 

(xO’yo)  5 Ray  intersection  with  x,  y-plane  , 

and 

9 = aspect  angle. 


where 


t^  = -sin  9 cos  a r + sin  9 sin  a <p+  cos  9 z 


cos  ot  = 


2 2 
xn+  yn 


The  low  aspect  angle  rays  are  terminated  when  the  ray  exits  from  the  plume 
or  the  value  of  the  exponential  damping  term  is  in  the  range 


10  ^10 


< -20  db  . 


(42) 


There  is  no  overdense  surface  ray  termination  of  cross  section  return  for  the 
low  aspect  angle  case.  The  cross  sections  are  not  range  gate  resolved  for 
low  aspect  angle,  but  they  are  Doppler  resolved. 
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3.  MAIN  PROGRAM 


A flow  chart  of  the  main  program  is  shown  in  Figure  4.  The  main  program 
performs  the  x and  z integrations  of  the  cross  sections  using  the  Trapezoidal 
rule  end  calls  subroutine  RAY  to  perform  the  integration  along  the  ray  (in  s). 

For  low  aspect  angles  (9  < $ ),  the  integration  in  z is  transformed  into  an 
integration  over  y (see  Figure  1)  by  a suitable  modification  of  the  step  sizes 
integration  limits. 

The  radar  properties  are  read  in  by  the  main  program , and  LAPPIN  is 
called  to  read  and  store  the  plume  properties  from  theAeroChem  LAPP  code.^ 

If  the  radar  aspect  angle  is  greater  than  90°  the  plume  is  reflected  about  a 
bisecting  x,  y-plane  (Figure  1)  and  the  aspect  angle  (0)  is  replaced  by  its 
complement  (180°  - 0). 

The  integration  step  sizes  are  set  equal  to  a fraction  of  the  larger  of  the 
radar  free  space  wavelength  or  the  jet  radius  by  an  input  data  card  (C  in  Figure  6). 

For  low  aspect  angles,  the  s and  z integration  step  sizes  are  interchanged. 

If  the  input  parameter  INCOH  is  set  equal  to  I.  the  step  sizes  are  all  set  equal 
to  the  step  size  in  z regardless  of  aspect  angle. 

The  limits  of  the  i.  integration  are  taken  to  be 


and 


ZOO 


PLRMAX 
tan  9 


ZUM 


PLRMAX  . ni 

uT7'  + PLZ(NZ> 


(43) 


(44) 
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Determine  Integration 
Step  Sizes  STEP,  DK  LX, 


Calculate  Cross  Sections 
for  Current  Step  in  X 
Integration 


Calculate  Z integration 
Limits  ZOO,  ZLIM 


.X  > XUM 


Print  Step  Sizes  and 
Integration  Limits 


Calculate  Cross  Sections 
for  Current  Step  in  Z 
Integration 


Overdense  Surface  Radii 
PLROVP  (I) 


ZO>TESTl 


Doppler  IUn  Maximum 
Velocity  and  No.  Bins 
PLVMAX,  NI)B 


Print  Cross  Sections  for 
Current  Range  Cell 


Range  Cell  I-ength 
FL  


7()>ZUM 


Print  Cross  Sections  for 
All  Range  Cells 


Figure  4 - Main  Program 


where  PLRMAX  is  the  maximum  plume  radius  and  PLZ(NZ)  is  the  length  of 
the  plume.  The  values  above  are  multiplied  by  tan  0 for  low  aspect  angle 
cases  to  give  the  corresponding  limits  for  the  y integration. 

The  plume  data  points  at  a given  axial  station  are  tested  for  the  overdense 
election  density  condition  (Subsection  2.5)  starting  with  the  maximum  radius 
and  stepping  Inwards.  The  radius  of  the  overdense  surface  at  the  axial  station  is 
taken  to  be  the  linear  interpolation  between  the  first  overdense  point  encountered 
and  the  preceding  point. 

The  number  of  Doppler  bins  required  for  storage  of  the  Doppler  resolved 
incoherent  cross  section  is  determined  by  searching  for  the  maximum  gas 
velocity  (PLVMAX)  projected  along  the  radar  line-of-sight  and  dividing  the  maxi- 
mum velocity  by  the  radar  velocity  resolution. 


NDB 


PLVMAX  * 
DOFVEL 


(45) 


The  velocity  resolution  (VD)  may  be  given  directly  by  the  radar  properties  read 
in , or,  if  the  velocity  resolution  is  read  as  0,  it  w»ll  be  calculated  from  the 
frequency  resolution  (A  fR) 

7 re 

VD  =17  4fD  • <46> 

The  x and  z integrations  are  performed  using  the  Trapezoidal  Rule, 

[»(«,)-.  (S)l  . (47, 

0 

The  cross  sections  are  printed  when  z is  equal  to  a multiple  of  the  range  cell 
plume  intersection  length  (FL).  The  incoherent  cross  sections  is  presented  in 
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Doppler  bins  with  a maximum  projected  gas  velocity  and  corresponding  frequency 

shift  printed  for  each  bin.  The  total  cross  sections  for  the  current  range  cell 

I 21 

are  printed  in  units  of  square  meters  (sm),  and  the  value  of  10  Log10  | <Vlm  | 
(dbsm)  is  also  presented.  The  total  cross  sections  and  the  Doppler  resolved 
incoherent  cross  section  are  summed  over  all  range  cells  and  printed  at  the 
conclusion  of  the  integration. 

Subroutine  Ray  is  called  to  perform  the  integration  along  the  s-axis 
(see  Figure  1).  The  arguments  include  the  ray  intersection  with  the  x,  z-plane 
(xQ,  z0).  The  intersection  points  form  the  grid  over  which  the  volume  integra- 
tions are  carried  out.  The  rays  are  formed  by  incrementing  the  intersection 
points  according  to 


and 


k - 1 
x 

^ dx 

k = 1. 

X 

2.  .. 

XLLM 

• ’ 4x 

(48) 

k - 1 
z 

) iz  * zoo  ■ 

k = 1. 
z 

2.  .. 

ZLIM 

* ’ Az 

(49) 

as  the  integration  proceeds. 

For  low  aspect  angle  cases  the  integration  over  z is  carried  out  with  limits 
and  step  sizes  appropriate  to  an  integration  over  y.  and  the  point  (x^.  zQ)  is  the 
intersection  of  the  ray  with  the  x.y-plane. 
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4.  SUBROUTINES 


4.1  Subroutines  LAPPIN  and  CYLDER 

Subroutine  LAPPIN  reads  in  the  plume  properties  from  thG  AeroChem 
LAPP  code/1)  The  LAPP  output  format  with  LAPP  FORTRAN  variables  are 
shown  on  the  next  page.  To  be  used  as  input  to  the  present  program,  the 
deck  must  be  preceded  by  a title  card  in  80  column  alphanumeric  format,  and 
followed  by  an  end  card  which  has  -1.  E30  punched  in  columns  75-80. 

The  plume  properties  are  given  in  the  x,  z-plane  (see  Figure  1)  at  the 
points 


(xj.z.)  = (PLR(I),  PLZ(J))  . 


The  plume  properties  stored  are: 


PPLE(I.J) 

PPLF(I.J) 

PPLV(I.J) 

PPLR(I.J) 

NZ 

NR  (I) 


_g 

electron  density  at  (x^z.)  (cm  ) 
electron  - neutral  collision  frequency  (sec-1) 
axial  gas  velocity  (m/sec) 
plume  radius 

number  of  axial  stations  along  z-axis 
number  of  radial  points  at  z. 


(50) 


All  distances  are  in  units  of  jet  radii  and  the  axial  gas  velocity  is  converted 
from  the  LAPP  units  of  ft/sec  to  Aerodyne  LAPPIN  m/sec. 

Various  consistency  chocks  are  made  on  the  input  data.  If  any  input  errors 
are  encountered,  the  variable  IER  is  set  equal  to  1 , and  a message  is  printed. 

The  values  of  electron  density,  collision  frequency  and  gas  velocity  are 
set  equal  to  1.  if  their  AeroChem  LAPP  value  is  less  than  1. 
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WBITT' 


CYLDER  can  be  cal?  3d  m place  of  LAPPIN  by  changing  the  parameter  ITEST  on 
a data  card  in  the  main  program.  The  test  plume  defined  by  CYLDER  is  a cylinder  of 
constant  electron  density.  CYLDER  supplies  a plume  data  point  every  tenth  of  a jet 
radius  in  the  radial  direction  and  integer  multiple  of  a jet  radius  in  the  axial  direction. 
The  plume  properties  are  set  by  data  cards  in  CYLDER  to: 

Length  = 2 Jet  radii 

Radius  = 1 Jet  radius 

8 3 

Ne  - 10  /cm  (electron  density) 
v - 1011  Hz  (collision  frequency) 

v =1.  (axial  gas  velocity) 

TABLE  1.  LAPP  CARD  OUTPUT 


(1.  El 0.3,  60X,  El 0.3) 


XORJ 

BIG 


1.  E30 


- axial  station  scaled  by  rj 

- card  type  1 identifier 


Card  Type  2 (1P8E10.  3) 


YOUT 

rik 

- radial  location  scaled  by  rj 

ECC 

Neik 

electron  density  (cm~3) 

XNEU 

‘'ik 

- electron-neutral  collision  frequency 

1 T 

Tik 

exhaust  temperature  (°K) 

1 PPCO 

PCO 

= partial  pressure  of  CO  (atm) 

PPC02 

pco2 

partial  pressure  of  C02  (atm) 

1 PPII20 

PIl20 

partial  pressure  of  H20  (atm) 

u 

uik 

exhaust  axial  velocity  (ft/sec) 

Molecular  species  partial  pressures  are  part  of  the  normal  LAPP  output  but  are 
not  used  in  the  RCS  calculation. 


Note:  Each  type  1 card  is  followed  by  a variable  number  of  type  2 cards  but  not  more 
than  25. 
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4.  2 Subroutine  RAY 


Ray  performs  all  of  the  calculations  necessary  for  a single  integration  ray. 

The  ray  is  parallel  to  the  radar  line-of-sight  (s-axis  in  Figure  1)  and  passes 
through  the  x,  z-plane  at  (xQ,  0,  zQ)  or  the  x,  y -plane  at  (xQ,  zQ,  0)  for  low  aspect 
angle  cases.  The  plume  properties  are  passed  to  RAY  through  COMMON/PLUME/ 
and  the  output  is  passed  through  COMMON/RAY/.  The  subroutines  called  are 
INTER.  SUND,  and  INDEX.  A flow  chart  of  RAY  is  shown  in  Figure  5. 

The  minimum  and  maximum  ray-plume  intersection  points  are  determined 
by  dividing  the  plume  into  cylindrical  slabs  as  shown  in  Figure  2.  A test  is 
made  at  each  axial  station  for  a ray-cylinder  edge  intersection.  The  y coordi- 
nate of  an  intersection  point  is  (from  Figure  2) 


y = -S  sin  B . (51) 

The  ray  is  parallel  to  the  y,  z-planc  so  that  the  radial  position  vector  R from 
the  z-axis  to  the  intersection  point  is  given  by : 


R = v y • 


(52) 


Therefore  the  S coordinate  of  the  intersection  [joint  is 

n/r2  - \ 


S - ± 


sin  B 


(53) 


where  the  sign  is  negative  for  SM1N  and  positive  for  SMAX. 


SM1N  or  SMAX  may  occur  at  a cylindrical  slab  face  (see  Figure  2)  rather 
than  a slab  edge,  therefore,  each  axial  station  is  tested  for  ray  intersection 
with  the  slab  face.  The  S coordinate  of  the  faee  intersection  at  z.  is 


( zi  ~ ) 

cos  B 


(54) 
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Each  intersection  S value  is  compared  with  previous  maximum  and  minimum 
values  to  determine  the  new  SMIN  and  SMAX. 

The  ray  coordinates  for  high  aspect  angle  (0  > 0Q)  are  given  by: 


S = SMIN  + (I  - 1)  STEP  (55) 

X = XO  (56) 

Y = -S  SIN  (0)  (57) 

R = \x2  + y2  . (58) 

For  low  aspect  angles  z and  R are  given  by: 

z = ZOLD  + TZ  STEP  (59) 

and 

R = ROLD  + TR  STEP  , (60) 


where  TZ  and  TR  are  components  of  the  unit  tangent  to  the  ray,  and  STEP  is 
the  integration  step  size. 

INTER  is  called  to  determine  the  plume  properties  at  the  point  (z,  R)  by 
interpolation  in  the  AeroChem  LAPP  plume  data.  INDEX  is  called  to  calculate 
the  index  of  refraction  using  the  interpolated  plume  properties.  SUND  calcu- 
lates the  integrands. 

The  tangent  vector  (t ) calculated  during  the  current  Integration  step  is  the 
tangent  of  the  ray  at  the  next  point.  Therefore,  the  currently  calculated  value 
determines  the  position  of  the  integration  point  following  the  next  point. 

Portions  of  the  plume  are  excluded  from  the  volume  integrations  for  the 
cross  sections  by  terminating  the  ray  according  to  tne  conditions  shown  in 
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Figure  5.  For  low  aspect  angles  (9  < 9^,  Subsection  2.  8)  the  ray  is  terminated 
when  the  attenuation  term 


exp  1-2  J ar.dS' 
\ SMIN 


is  less  than  -20  db  or  when  the  value  of  IOUT  is  equal  to  1.  Subroutine 
INTER  returns  a value  of  1 for  IOUT  when  an  attempt  is  made  to  extrapolate 
the  plume  data  to  a point  outside  the  plume. 

High  aspect  angle  rays  are  terminated  by  IOUT  = 1 or  when  the  ray 
encounters  the  overdense  surface  defined  in  Subsection  2.5. 

The  Simpson's  rule  integration  scheme  for  the  ray  integration  of  the 
cross  sections, 


J f<x>  dx  = (f2  f 4fl  + fo)  • 


requires  that  the  ray  terminate  on  an  odd  number  of  integration  steps. 
Furthermore,  the  attenuation  and  phase  portions  of  the  integrands  require  the 
evaluation  of  integrals  at  Simpson's  rule  midpoints  (even  numbers  of  integration 
steps).  This  is  accomplished  by  using  the  "3/8  rule," 


a 

/ {(x>dx  " idx  (f0+  3fl  + 3f2+  f3  ) 


to  evaluate  the  integrals  in  the  attenuation  and  phase  (P  in  Subsection  2.4)  for  even 

number  integration  steps  and  Simpson's  "1/3  rule"  for  odd  number  steps.  This 

5 

ensures  tliat  the  combined  integration  has  a truncation  error  of  dx  . 
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The  incoherent  cross  section  contribution  for  each  integration  step  is 
stored  in  a "Doppler  bin"  dependent  on  the  local  axial  gas  velocity.  The 
Doppler  resolved  contributions  from  the  incoherent  cross  section  are  deter- 
mined by  a Trapezoidal  integration  over  each  integration  step. 

The  coherent  and  overdense  cross  sections  for  the  current  ray  are  set 
equal  to  zero  if  NFACE  equals  1.  The  first  section  of  RAY  that  determines 
SM1N  and  SMAX  will  set  NFACE  equal  to  1 if  the  ray  intersects  the  plume- 
rocket  interface  at  z = 0. 

4.  3 Subroutines  INTER  and  PROF 

Subroutines  INTER  and  PROF  perform  interpolation  of  the  plume  data  to 
obtain  the  values  of  Ne,  v , and  v at  the  point  R,  Z.  A linear  interpolation  in 
Z is  performed  by  INTER  and  a quadratic  interpolation  in  R is  performed  by 
PROF. 

INTER  finds  a pair  of  axial  stations  such  that 


< z < z 


i+1  * 


PROF  is  called  for  each  axial  station  to  perform  a radial  interpolation  obtaining 
the  plume  parameter  values  at  (R,  Z^)  and  (R,  j ).  INTER  performs  a linear 
interpolation  in  Z between  the  axial  station  values  to  give  the  plume  values  at 
(R,  Z).  Points  outside  the  plume  generate  exponentially  decaying  plume  values. 
The  attenuation  factor  is  the  inverse  of  the  free  space  wavelength. 

4.4  Subroutines  INDEX  and  SUND 

INDEX  determines  the  complex  index  of  refraction  (N  - XIND  - i Y1ND) 
from. 

Re  \Zn2  = P ♦ yy2  + xx  ) 
and 

ta  = -J~  (Vxx*  ♦ yy2  - xxj 


(64) 

(65) 
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where 


Ms)1 


Ii  INDEX  = 0 the  index  of  refraction  is  set  to  the  free  space  value  N = 1. 

SUND  calculates  the  integrands  from  the  interpolated  values  of  N ,v  and 

9 - 2 e 

the  index  of  refraction.  The  parameters  p“=]/ 2(1+  r ) and  cos  fi  are  set 
to  1 and  -1  respectively.  The  output  from  SUND  is: 


(Subsection  2.3) 

a (attenuation  coefficient,  Subsection  2.3) 
r 9*1 -1/2 


DTAU 


= dr  --  reNc  _1  f (~) 


2N  (j/c  (Phase  Integrand,  Subsection  2.4) 


5.  INPUT/OUTPUT  FORMATS 


5. 1 Input  Format 

The  card  deck  for  a run  using  a FORTRAN  source  deck  on  a C DC  6600  with  a 
SCOPE  operating  system  is  shown  in  Figure  6.  Card  B gives  the  radar  properties, 
Card  C gives  the  Integration  step  sizes  and  Card  D sets  the  option  switches.  The 
AeroChem  LAPP  data  defining  the  plume  is  described  in  Section  F.  The  required  plume 
data  card  format  is  shown  in  Subsection  4. 1.  Card  G signals  the  end  of  the  plume  data 
cards.  Cards  A-D  may  be  repeated  in  Section  H for  additional  runs  using  the  same 
plume  data 

CARD  A:  RADAR  Properties  Title  Card  (20A4) 

CARDB:  RADAR  Properties  (Free  Format) 


9 (deg) 

Aspect  angle  (0°  for  nose  on  viewing) 

f (Hz) 

Radar  frequency 

i-j  (m) 

Jet  radius  at  matched  pressure 

Range  (m) 

Missile  range 

RESOLUTION  (m) 

Radar  range  cell  resolution 

Divergence  (deg) 

Radar  beam  divergence 

Doppler  Velocity  (m/sec) 

Radar  doppler  velocity  resolution 

noppler  Frequency  (Hz) 

Radar  doppler  frequency 

Resolution  (set  to  0.  if  velocity  resolution  is 

given) 

Low  Aspect  Angle  (deg) 

Angles  less  than  or  equal  to  this  angle  are  low 
aspect  angles  treated  as  described  in 
Subsection  2.  8 

Gamma 

Mean  fluctuation  factor 
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CARD  C:  Integration  Step  Sizes  (Free  Format) 

FRACTS  Fraction  of  free  space  radar  wavelength  chosen  for 

integration  step  size  along  the  ray 

FRACTX  Fraction  of  free  space  wavelength  chosen  for 

integration  step  size  along  x-axis 

FRACTZ  Wavelength  fraction  used  for  integration  step  size 

along  z-axis  for  high  aspect  angle  cases  and  y-axis 
for  low  aspect  angle  cases 


CARD  1):  OPTION  Switches  (Free  Format) 

ITEST  0 Plume  data  from  card  deck  in  Section  6,  Figure  6 

1 Constant  electron  density,  cylindrical  plume  data 
from  subroutine  CYLDER 

IPRINT  0 No  printout  of  integrands  every  Nth  step  in 

ray  integration 

N Printout  of  integrands  every  Nth  step  in  ray 
integration 

IATT  1 Integrand  attenuation  factor  attenuates  cross  sections 


0 No  attenuation,  attenuation  factor  set  equal  to  1 

ID  EX  1 Index  of  refraction  calculated  from  local  electron 

density  and  collision  frequency 

0 Index  of  refraction  set  equal  to  1 


lOV!) 


Uav  integration  terminated  when  overdense  surface 
encountered  for  high  {$  > 15  ) aspect  angle 

No  overdense  surface  ray  termination 


0 


REPEAT  CARDS  A - D FOR 
ADDITIONAL  RUNS 


SYSTEM  A .47  KI'T 

0,0.1. 1.1.0 


1 --> . . 1 . :*2  E9 . . 4 . 5 . K4 . 1 5 . . . «5.  .400 . , 0 . . 1 r> . , 1 . 
lf>  DECREES 

* EOll 


FORTRAN  SOURCE  DECK 

* COR 

LOO. 


FTN. 

^ KICKR,  T2000,  CM  120.000.  227s 

Figure  6 - Card  Deck 

.10 

I 

1 

I 


tait 

AL-1406 


IIIUFF  1 Surface  roughness  factor  and  applied  to  overdense 

cross  section  return  (Subsection  2.  5) 

0 Surface  roughness  factor  set  equal  to  1. 

CARD  E:  PLUME  DATA  TITLF  CARD  (20A4) 

SECTION  F:  Refer  to  Subsection  I.  1 - Table  1 . For  a constant  cylinder  test  case 
with  ITEST  1 on  Card  D,  Section  F is  omitted. 

CARD G:  END  OF  PLUME  DATA  0.  in  columns  1 and  2,  -1.  E3  in 
columns  7'i— 30 

SECTION  11:  Repeat  Cards  A-l)  for  additional  runs  using  same  plume  data 
in  Section  F. 

5.  2 Output  Format 

Appendix  B shows  the  PARCS  output  lor  Svst(  A using  the  input  data  of 
Figure  (i,  Subsection  f>.  1.  The  first  page  of  the  output  shows  the  plume  data  used  in 
tlie  output  subheading  "SYSTEM  A :!7KFT"  taken  from  plume  title  card  E (Figure  6). 
'Hie  data  from  input  data  cards  \ through  I)  is  printed  below  the  plume  title  cart! 
"SYSTEM  A :!7  EFT"  in  the  same  units  as  those  shown  in  Subsection  5. 1.  The  last 
line  of  the  first  output  page  shows  the  tower  and  upper  limits  of  the  y integration 
ZOO  and  ZLIM  in  units  of  jet  radii  DU  FT).  Tin*  low  aspect  angle  definition  0^  and 
tlie  radar  aspect  angle  o are  !>oth  equal  to  15  ('eg  so  Lhat  tills  is  a lew  aspect  angle 
ease  (Subsection  •>.  »).  Therefore,  die  integration  in  /.  (along  the  plume  axis. 

Figure  1 ) is  transformed  into  an  integration  along  the  v-axis.  The  number  of  integra- 
tion steps  taken  along  Ute  y-.uxis  is  given  by  NZM  and  DZO  is  the  y integration  step 
size  in  units  of  jet  radii  (RJKT).  If  this  had  been  a high  aspect  angle  case  (0  > 0q), 
YAK),  ZLIM,  NZM,  and  DZo  would  a!)  have  referred  to  tlie  /.-axis  and  there  would 
have  been  no  integration  along  the  y-axis.  DXO  is  tlie  step  size  along  the  x-axls 
and  STEP  is  die  integration  step  size  along  the  my  (Figure  1).  Both  DXO  and  STEP 
are  in  units  of  Jet  radii. 
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The  second  and  third  pages  of  the  sample  output  show  the  cross  sections  for 
that  portion  of  the  plume  in  the  first  range  gate  (Subsection  2.  6).  The  range  cell 
number  and  length  are  shown  on  the  first  line  of  output  page  B-2  as  1 and  10*®  jet 
radii.  The  plume  is  contained  entirely  within  the  first  range  cell,  but  if  additional 
range  cells  were  required  to  cover  the  entire  plume  their  printouts  would  follow  one 
another  in  order  after  the  first  range  cell.  The  incoherent  cross  section  of  each 
range  cell  is  "Doppler  resolved"  (Subsection  2.7).  The  first  "Doppler 
bin"  (page  B-2)  gives  the  total  incoherent  cross  section  for  those  portions  of  the 
plume  where  the  plume  axial  gas  velocity  projected  along  the  radar  line-of-sight 
(ray.  Figure  1)  is  less  than  or  equal  to  the  "DOPPLER  VELOCITY"  of  300  m/sec. 

The  second  Doppler  bin  contains  contributions  from  regions  of  the  plume  where  the 
projected  gas  velocity  is  in  the  range  of  300  to  800  m/sec.  The  "FREQUENCY  SHIFT" 
is  the  corresponding  Doppler  frequency  shift  of  the  radar  wave  in  hertz. 

The  cross  sections  are  shown  in  two  systems  of  units;  "sm"  or  square  meters 
and  "dbsm"  where 


(r(dbsm) 


10  Log 


10 


a (sm) 

1 m2 


(72) 


The  total  cross  sections  for  the  first  range  cell  are  shown  on  page  B-3.  The  total 
incoherent  cross  section  (SINC)  is  the  sum  of  tire  contribution  from  each  Doppler 
bin.  SCOH  is  the  coherent  cross  section,  SOD  is  the  contribution  from  the  overdense 

O 

surface,  and  RCS  is  the  total  radar  cross  section  given  by  the  sum  (in  units  of  m“) 
of  the  incoherent,  coherent,  and  overdense  cross  sections.  Output  pages  similar  to 
pages  B-2  and  B-3  will  appear  for  each  range  cell. 

Pages  D— 1 and  B-5  show  the  total  cross  sections  for  all  range  cells.  Hie 
Doppler  bin  incohe—mt  cross  sections  are  equal  to  the  sum  (in  m2)  of  the  cross  sec- 
tions from  the  corresponding  Doppler  bins  in  each  range  cell.  The  total  cross  sec- 
tions for  all  range  cells  and  Doppler  bins  are  shown  on  page  B-5. 
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The  "Day  File"  on  the  last  page  shows  a CP  (central  processor)  time  of 
"1850. 4b  see"  f >r  the  CDC  GCOO . phis  run  would  Lake 

:>.4  x 1x50  = 6,2!H>  sec 
on  a C1K  G 1 Of » computer. 

5 . ! Program  l so 

In  order  to  aseer'cin  the  accuracy  of  the  PARCS  calculations,  a number  of  test 
eases  (for  which  analytical  solutions  are  available)  are  employed.  These  test  cases 
also  permit  a rapid  ansessnu  at  of  the  effects  of  varying  stew  sizes  employed  in 
integrations  along  the  rays,  As  and  eh  uiging  the  total  number  of  rays,  i.  e. , changing 
Ax  ami  A/..  The  general  eonelusion  < these  tests  cases  is  that  the  incoherent,  under- 
done • RCS  is  least  sensitjv  to  step  size  c hanges.  The  incoherent  RCS  requires  five 
steps  per  wavelength  to  give  a '!  do  ag>vemenl  with  the  analy deal  solution.  At  least 
10  steps  per  wavelength  are  nculed  for  a 2 db  agree.ment  uitli  the  closed  solution  for 
the  coherent  overdun.-e  and  und<  rleum-  rcRmis. 

Figure  7 -hi iv  s Uu  > ckwla'co  1 -band  underdense,  incoherent  RCS  > . a function 
of  aspect  angle  for  a uniform  cvlinoer  0.  xm  long,  0.  Im  radius  (R  ) having  an  electron 
lens  t_\  -.i  l'i*  ,-em  ’ and  'hsi  i ' •.  • aetu  v of  10  ^/see.  There  is  exact  agreement 
at  00  with  the  conation  for  the  eoilision  corrected  first  order  Horn  approximation 

to  the  Incoherent  'moss  .vui'"!* 


a — y N“  <f»V  , (73) 

1 * i t - r 

> u 

where  V is  the  yoluv.e  and  the  other  terms  are  given  in  Section  2.  At  lower 
as f>eet  angle-,  here  is  a decrease  in  tin-  RCX  due  to  increasing  path  lengths  and 
resultant  attenuation  in  the  plasma. 

h 1.  farcin,  ' Radar  Reflevtivttv  of  Turbulent  Rocket  Kxhaust  Plumes  — Derivation 
of  Kquauo-is, " A'-r'Oi<  ;u  Research  ! alxjralories,  TP-167  Supplement, 

Nm  i mlicr  ItH  ’.i,  i*.  yo 
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Figure  7 - Underdense,  Incoherent  RCS  for  Uniform  Cylinder  (ne  = 10*0 
V - lOH  sec-*  L = 0.  8m,  R = 0.  4m) 


The  underdense,  coherent  return  for  the  same  cylinder  is  shown  in  Figure  8. 
The  closed  form  solution  is  obtained  from'  ' 


a - 


32rr2r2L2R4N2 
e c e 


sin  (kL  cos  B ) 

2 

sin  ( 2kRc  sin  B)~  j 

kL  cos  $ 

(2kRc  sin  $ )3^2 

(74) 


At  90  there  is  again  good  agreement  with  the  exact  solution  for  all  steps  sizes. 
However,  at  aspect  angles  less  than  about  70° , it  becomes  necessary  to  take  at  least 
ten  steps  per  wavelength  to  attain  less  than  5 clb  error. 


The  PARCS  calculated  return  from  an  overdense  smooth  cylinder  of  finite 
conductivity  is  shown  in  Figure  9.  The  closed  form  solution  (solid  line)  is  calculated 
from^ 


with 


J (N-  k)  eak  ‘ r dS 

S 


y/  N2  + N2 


(75) 


(76) 


where  P2  is  the  reflectivity  (Subsection  2.  5)  and  the  surface  integral  is  evaluated  on 
the  surface  of  the  cylinder.  It  is  necessary  to  employ  greater  than  10  steps  per 
wavelength  for  the  overdense  cross  section  in  order  to  get  good  agreement  with  the 
closed  solution. 


5.  J.  Jurem,  "Radar  Reflectivity  of  Turbulent  Rocket  Exhaust  Plumes  - Derivation 
of  Equations,"  AeroChem  Research  Laboratories,  TP-167  Supplement, 
November  1969,  p.  50. 

6.  Methods  of  Radar  Cross  Section  Analysis,  ed.  J.W.  Crispin  and  K.M.  Siegel; 
p.  105,  Academic  Press.  1968. 


35 


37 


Figure  9 - Test  Case  for  AFRPL  Overdense  RCS  Routine 


The  total  coherent  and  Incoherent  cross  sections  for  SYSTEM  D at  23.9  km  was 


calculated  for  various  step  sizes  using  the  PARCS  code  and  the  input  data  of  Appen- 
dix C.  The  results  are  shown  below. 

Step  Size  Total  Incoherent  Total  Coherent 

FRACTX.FRACTS,  Cross  Section  Cross  Section  CPU  Time  (sec) 
FRACTZ  fdbsm)  (dbsm)  CPC  6600 


0.1 

0.1 

0.5 

-45.5 

-59.7 

100 

0.2 

0.2 

0.2 

-45.9 

-69.9 

84 

0.1 

0.1 

0.1 

-45.9 

-67.7 

490 

0.05 

0.05 

0.05 

-45.9 

-67.7 

3,212 

All  step  sizes  must  be  taken  as  one  tenth  of  a radar  wavelength  to  get  convergence 
for  the  coherent  cross  section,  but  the  incoherent  cross  section  is  close  to  convergence 
with  FRACTZ  = FRACTX  = FRACTS  = 0.  2. 

The  CPU  time  is  approximately  a linear  function  of  the  step  size  in  each  direction. 
For  example,  decreasing  the  step  size  in  the  z direction  by  a factor  of  five  in  going  from 
the  first  to  the  third  case  in  the  table  above  increased  the  run  time  by  a factor  of  4. 9. 
The  CPU  time  on  a C DC  6400  computer  would  be  a factor  of  3.4  times  longer  than  for 
the  CDC  6600  (see  Appendix  C,  System  D). 
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APPENDIX  A 


PROGRAM  LISTING 


.-G_ 

r 

r. 


flFPODVM c RCS  PROGRAM  CY=2C1 

PR'GPAM  P APTS  (OUTPUT,  TAP "6=  OUTPUT, TA®E1 , TA°E3 ) 
C0NH0N/PttiME/NZ,F4.2«tfCHNRTluv»  ,PL«mfr2,2Si  ,PL£(lfru ,25) 
*,PLP(1jO , 25) , I ATT , IDEX, IOVO,PLV (10„ ,25)  ,NNR(13C),  IRUFF, 

V?5*-  rPPtr^rt  « Jr-^f 


COMMON/RA  v/STUOA  , TTILOA  (?)  ,f MUOQ(2)  ,DOPSTE(5? ) ,PlROVD(1CO  ) ,GA MM  A 
CPHMOH/ AN€L £/ THE T A, T«tT AO, T Hr T AO, TH5 TOO, COS TH-  , SI  NT Hr  ,TTME TA 
DIMENSION  TITLE(ZC) ,TITLEi(2j) 

OTMFNSION  TGOH<2>  , 5MUOO 42  > « XCOH < ?T  rX£HH?f  , 

1 OLOCHX(2)  ,OLOOOX  (?)  ,OLOCH7  (2)  , OLDODZm  , 

» COai-NO (5  (?)  ,00PX(5.;  ) , CL  ODOR  (5  [• ) ,01-00^?  f 9-? ) ,005000(50)  , 


r 

r 

-C- 


T°OP  (c0  ) ,D!’TO OP  ( 5C  ) 
pcai_  (.AMROA 

nATA  nrG/.C174E329?5/,DI/3.  1-*1S9?6j4/,C/2. 393E  8/ 

IP? INT  = N FOP  INTEGRANDS  PRINTED  EVE® Y NTH  STEP  IN  S 
4-A.z T— - r-  4—c Oft  ■ CALOULA  T F 0 A-Y  TE-NOA  PEON 


-C_ 

r 

r, 

r 

c 

r 

-C- 

r 

r 


? CCR  NO  ATTENUATION 

IOpx  = 1 PC-O  CALCULATED  INO£X  OF  c-FRACTION 
D r0R  INDEX  = 1 

TCV-9  = 1 POO  OVER-DENSE  SURFACE  RAY  TPRMINATIGN 
1 N'>  OVEROENSE  CUP  FACE  ° AY  T ERMI  NAT  I ON 

aVRaO£NSF~SCFFACr  ROUGHNESS  G AL-GU-L  AT  ION 

? NO  SURPACE  POUGHNES1'  CALCULATION 
I TP  ST  - 1 CYLINDRICAL  PLON*  OF  CONSTANT  DENSITY  PROM  CYLOE® 
r PL UME  DATA  EP OH  LAPPIN 

PO40TS,POAr,TX,FPACTZ  * FRACTION-  OF  FRE-  SPACE  MAV£LEMGTH  CHOSF  N 
<■  )R  INTEGRATION  STEP  SIZT 

TA-O-iV  ^4, AaOPST  - ASPECT  ANGLE  -USING-  LOW-  ASPECT  - ANGLE1  -SCHENE  MuES^FFST 
PUNCTinn  state Mf NTS 


FOn<XX,YV,Z?TOM>  s XX  - ZZMW  4-  Y Y * Y Y / (CM*  ON ) ) X ( 1 • ♦ ,S*YYVOu) 

D9rTN(X) = i9,#ALOG13 (AHAX1(X,1,E-1P)  , 


NTTH"'-  = C 
1P-C-  CONTTNUP 

r> 

C RP AO  AND  POINT  PAOAR  OAT A 

r* 

8CA0  <Gt300)  TITLE 

IP  (PCP(5*,NE,r.)  STOP 

RFAO(t,* ) THPTA0»BF0?O,«Uc»  »RANG£  ,R€SOLU,OiV£RG, 

1 prpvrL,D0PFOQ,THET:r  , GAMMA 
OF  AO  < 5,  • » FRAG  T S»  FR  AG*X  .FRAGT? 

AO <*(,•)  ITEST,IPOINT,IATT  ,IpPX,TOVO,  IRUFF 


A-I 


r pr an  jm  PLUM-  DATA 


Ir I NTIMES .NF. 0 ) GO  TO  30 

IEHME^!)  19,29,44— — — 

1C  CALL  PYLOERdFR, TITLED 

GO  T&—t6 — 

2*1  CALL  LAPPIN(I£P,T ITLE1) 

— 3£  IF4IfR.NE„P4  STOP- 

NI-IH£S_=_ 1 

WPTTF  (F,  205) 

__  295  EO-CHA T ( »4-AER00*N£  J*GS  -PP-OGP-AH* /1‘  - 

WRITE  *•-,  20*  > TITLE1,  TITLE 

...iO:  FGFMA-TI2CAA4- 

W»ITF  (£,  2 20  ) THETAO,RFREQ,RJET,RANGE,RESOLU,OI «ERG, 

1-0-OP-4E.U-OO  EE.P.Q  , THE  LC-Q^-GAMMA 

9 2r.  FO°MAT(*0*»T5»*THETAO*»T15,*RFC?EQ*»T25»*RJET*,TJ5  , * R ANGE* , T hS , 

1 *R -SOUJT I GN*,T  57, *GI MERGENCE*  / 5 X , 1 PGG1 9,3/ »J-»,T5,  * OOPVEL*-,  T 1 

:*OOPFOQ*,T?c,  *THETOO*,T35,*GAMMA*/5X,lPl»Gia.  3) 

— — MR  I TE  t5^222-T  4-PACTG,FPACT  X,FRACT7  

222  FO-MAT(*0*,T5,»FRACTS*,T15,  *FPACTX*  ,T25,*F=?ACTZ*/5X,  1P3G10  . 3) 

WUl£-lL+Z2kl-~ I-I.£S-I.,4P-R.IM.T-^I-AJ-T.,.Xa£)C.,.X.OV.Oy4R-Ug-g 

224  FOPHATC*0*,5X,*ITFST*,T15,*IPRINT*,T25,*IATT*,Tl5 ,*I0=X*,T,5, 

1 *IOVO*,I55,*IRUFF*/feI44)  - 

r 

r.  PLUME  REFLECTED  ASOUT  -BISECTING  XV  PLANE  FOR  THETA  GT  9j  DEGREES  - 

r 

DO.  of..  Tsuaz 

NMAP  - I 

IF  ( THEIAD.GI*  9w)  NMAP..^_N7  . 1-- *-4-  

NR  T = NNR(I) 

DO  50  J=1,NRI 
PL  (NMAP, J)  = PPL  F ( I , J) 

PLF INMAP, J)  s PPL F 1 1 , J ) 

PL  y/  (NMAP,  J)  = PPU/(I,J> 

PL-R  INMAP , J)  =.  PPLR  C I,J4  - — 

NR(NMAP)  = NNR(I) 

5Z  CONTINUE 
F ? CONTTNUP 

TF  ITHFIAn.GEdR^l  I HF  Tan  z:  mr.-TH-TAD 

r 

C MAXIMUM  PLUME  RADIUS  PLRMAX 
r 

PLRMAX  = C. 

OO  65  1*1, N7 

PI  SMA.X  ■*,  AMAJLj-LELL RKAJUPl  R I,  I ,UE  1 1 DJ 

f>c  CONTINUP 


A-2 


C AREA  - PROJECTED  AREA 

C OMf.G  = SAOJLR  FRE-OUENCX-  IN  -RAO IANS /SEC  — - - - 

f'  th^TA  = ASPECT  ANGLE  IN  RAOIANS 

r.  DKG  1 GTEP-JU-S  INTEGRATION  STEP  .IZES-t-N  X ANO  S IN  UW-S-Qf  RJ£A 

P n?:  is  THE  STEP  SIZE  IN  THE  AXIAL  DIRECTION 
C . F ° A C S = TRACTS 

FRAC7  = FRACT7 
FPACX  - FP  ACT  X 

IF < JjiE.IAO.-GT.  T-HE-IX-Q.)-  GO  TO-  70 

FOACS  - FPACTZ 

FRAC7  - TRACT  X - - 

T DMT,  - RPRED*2.*PI 

LAMiJOA  = 2.*PI*C/0N£G/RJ£T 
STro  = LAHPOA*FPACS 

OX-.-— —LA  NgQA*  F-RAC-X 

DZ*  - LAMBDA*FRACZ 

if (lamdoa-i.t  - - — 

fro-  FT--P-  - -P-*ACS*RJET 

ox:  = PRACX’PJFT 

OZ"  - FPAC7-»PJET - - 

an  TH'TA  - THETAO*OEG 

T H c T A C-  = TMET9IWT€f,”  * 

T (THETAD.LE.89.)  GO  TO  22E 

T-Ti»c^  »■■-*■  -iT-ec-M-e 

GD  TO  228 

22E  THETA  =■  TANTTHETA) 

2?"  C 0 S T H c = COS(THETA) 

S TNTwE  * SORT-H.-COSTj+E*CC^THE) 

AS-AOO  = 2. ♦PXO*OZO*R J"T**2 

-A-P-rA—  » apc>00«5INTHE 

IF  (THETA. LE.THFTAO)  ARPA  = ARE  A OD 

FF r*0  = *.  !*.F-1P*0MEG*0“EG*SINTHE*SINTHF 

RAMGE  = RANGE/RJET 

RPFOLU  * RPSOLU/RJET 

0 T 7PPG  = C I VE RG*D  EG 


r 7LTm  = ypprp  limit  OF  7 INTEGRATION  CHOSEN  TO  INCLUT;  ;.N0  SECTION 

p N7m  - humo-p  OF  POINTS  IN  AXIAL  DIRECTION"  — 

r 

7 0'  * . 

Dr '-7  = 0. 

7L-IH  ■ PLZ(NZ) 

Ir (THPTAD  .GE, 89. ) GO  TP  2L2 
IF(TjiRT*o,Lf  ,i,j  -GO-  TO 
Pr L 7 = PLPMflx/TTHET A 

79  '-  * -1 . *OcLZ  

GO  m 2U9 


A -3 


zt  ht- *-tHtRha*- - 

7 D : = -l.*PLRMAX 
GO  TO  - 

?tn  ZLIM  = 7LIH  ♦ DEL Z 

- TP  < THP  TAtt-GT -tT-H-E  TA-6F  -GO-  -TO-EGG 
ZLIH  = ZLIH*TTHETA 
70?  » ?T?0*TTH;T  A 


?i,s  N?m  = TNT  (1.G001M7LIH-ZG0) /OZG)  ♦ 1 

2E3  FORMAT  l*C*,T5,*Z0C*.T15**ZLI:l*,T2E,*NZM*f  T 55, 

i *ftKf*rT»5TJtPZfr*f T55'i*ST«»*y*- -*rlf^tfrrft-Kr»»T»t7r 
IF  (IPPINT  ,GT.  31  WPITl(E,25EI 

j c,  s-  -p T P I-»*  /* f>*  t T C i * ‘j  t R J«  f T g 0 » ♦ EX  » , 

1 T-5,*FTHnA*tT50  ,*TTILQA  (1 I * , T63 ,* TTIL OA < 2» * . TdC  , 

2 *FMUOOfl-H*,T^«*fMUOef?t-*v 


T-HT5  tn0*>  OETEPMTEtES  T+fE  -RMH-OG  -Of 
AT  r ACH  AXIAL  STATION 4Z  = PLZCTIt 


T-+F 


I = AXIAL  STATION 

00  3C " 1=1, N7 

-NR4-*  NR-f-T-F 

00  29S  JJ  = 1 »NRI 
0-—  NftI GO  »■■■  1 


J = ° AOIAL  STATION 


r 

A 


TfST  POO  OATH  At  STWtTI-ew-lNOPeH^-NTt^-^fTOW  OVf R0ENS5  SttREACS 


iF(j.n.i»  r.o  to  250 

I P < F (WW  € < I , M ♦ 1 , Jf  ,-Ff  OG.ONEG » >-  2 RE  r i 

zy  continue 


LTNPAt>  INT-PPCLATICN  BETWEEN  CURRENT  ANO  PREVIOUS  RADIAL  STATION 
OPIP-aMlNCP  PL  feftW 444 


PO-'TST  « FOOTPLET  I , J»  ,PLf  TI  ,J»  ,PFOO,OMt(i»  - 
1 r"f)<OLF  <!,  J-l)  ,PLFa,J-l),PPOn,OMFG» 

IF {FOOTS!, LT. 3.9" 14  GO  T-0  -? 5v 

PL^OVOd)  = FOOfPLE  <1  ,1  » ,PLF(I  , 1)  ,FFOO,OHS',»  • 

1 <atai.T-,  j-u  - » 

p ( rC0  <nLr  < I » J ) , PL  F < I , J > ,FF00*0MEG)  - 
X FnoiPLP<IW-l>  ,»LF{  I,  J-|  4,  Ff  00,  ONFGM 
GO  TO  TP3 

?«r,  Ptpovnut  = plp<i,j> 

GO  TO  *01 


A-4 


T"?  CONTINUE 


r OnPOLE^  YFLOCITY  RESOLUTION  1 MO.  DOPPLER  EINS 


NT  P - . 

IF  < TUFT  AO. GE. SR.  > CO  TO  IU< 

IF (OOcVrL . EQ. C . I OOPVEL  = DO0 F° Q*PT *C/ OHEG 

PLVMAX  - r . — - — — 

DO  31?  1=1, NZ 

NR4  - NR-t  I » ------  - - - 

DO  31?  J = 1 » NP I 

PL)/ MAX  = AMAXl<PUVMAX,RLV(I,jn 

31n  CONTINUE 

31?  CONTINUE  — •_  

oLv/MAX  = PL  VMAX4COSTHF 

NO”  = TNT  Cl, 0301* PLVMAX /GOP  VE-L4-+  1 - 

IP(N0n.LE.5CI  GO  TO  Ziu 

NaXX£Xe^-IA&4 

■»lr  PO°MAT  TOO  MANY  DOPPLER  CINS  FOR  ARRAY  STORAGE  ALLOCATE*) 

r> 

0 CALCULATF  IMF  RANGE  CELL  Lf NOTH -EL  

r 

u-  a M-IO-I-M  - - -ftANGi  IQI-VE-RC 

TANGLE  = ATAN(9W10TH/PrS0L'J) 

IF/TTH-- TA.LF.TANGLEl  GO  TO  3ZG 
FL  = «H'IDTH/SINTHC- 

GO  TO  -.315 

pl  = P'SOLU/COSTHE 

..  ^,T3!^-^E-CJMRTA-.Lg.TMETAG»-  EL-  = !.€■«* » 

r 

p THE  ")C  L0OOS  RANGE  OVER  NZ*  Z-VALUf  SUXIAL  STATIONS/  ANO 
p N»-  X-VALUrS (PA DIAL  STATIONS!.  THE  TRAPEZOIDAL  RULc  IS  USEO 
P FOR  T«r  x ANO  7 INTEGRATIONS, 

r 

GTMC  - 6. 

Tcr  OH  = 1. 

tsoo  = o. 

TOTNC  » r.  , 

TSTG  = 0. 

Trnmij  - n, 

TPPH(2)  = 3. 

smjoom  - t». 

S**lJ00C?>  = c. 

- no  »l|HOq 

OOPJNP i I ) = 0. 

tow»< n ^ c . 


A-5 


332  CONTINUE 

KK — * i 

KKK  = 1 



T'_cTL  - FL  ♦ ZOO 

r * 

00  <*?:  KZ=1,NZH 

r **** 

KC HUNT  = KCOUNT  4 l 

xing-«- 

XCnH(l>  = o. 

XOOH(?I  « •?. 

xnrdt  = c. 

xon<?>  =.  r. 

DO  334  I = 1 » NO  P 

nQvXri  u 

334  CONTINUE 


ZC  = ZCO(K7-l»  *070 


r CALCULATE  tHE  v/ALUFS  OF  XLIH  AM  0 NX  M 


34"  JF(KK.CT.N7)  GO  TO  3SC 

IF<ZC-,LT,PLZRK*Of  00  TO  35'- 
XL^m  = °L  P( KK , NR ( KK  ) ) 

Nxw  - TW*T4^0?1**4^H/3X44  -4~i— 

KK  = KK  ♦ i 

lEtKK.f.T.H?)  GO  TO  W 

no  to  74-) 

3 no  TM7<:  .L^PtZdrW  ♦***  « 

IF  (TH'TA.LE.THFTAO)  NXH=INT  (1.(J001*PLPHAX/OX:)  ♦ l 

DO  3*0  <X=1,NXN 

JZ uu 

X’  - (<X-1)*OXO 


r °»Y  CALCULATED  the  integrals  in  thf  s direction  ALONG  the  RAY 


CALL  PAY( XI  ,Z0  t OHEG ,RJ=  T,5T  EP , I PR  IN T , OOPVfL , N0» . L AMBO A» 


IF»<x.nT.?»  GO  To  373 

OL-OTCX  * 3. 

OL^CHXdl  = 0. 

OtnO*OU£4—»  

OLOOOXUI  = 0. 

OL  nOOX  <?4  ■ 0. 

00  36?  1 = 1, NOP 
OL  OOOPM  4 -*  


7?n  XINC  = XINC  *■  .5*  (STILOfl  ♦ OLO ICX > 

XRUHd)  -=  -XGDH4L)  .&MTTH.-0A-M!  ♦ 0L£€*XI444 

XC0H(2)  = XCOH (2)  ♦ . ?* ( TTI LO  A ( 2)  ♦ OLDCHX(2)) 

— xooo  - xaocii  * .s-mpmuodu > ->  oldogxtiu  — 
XOn<2>  - xn0(2»  4-  . 5*  (FMUOQ  (2 ) ♦ OLOOOX(2)) 

no  372  1=1, NOB 

nonx(D  - nopxd)  ♦ .5*  (oopstr<  i>  ♦ olooop(I)) 

OLOGOPII)  = DOPSTRU4 

7 7'’  CONTINUE 

GtniCX  ^ s T-ILOA  * 

OLnCHX(l)  - TTTLDA(l) 

QlflCax.<21...=  T TILDA <24 

oLocnxm  = pmuoo(I) 

0LP0DXI2)  = F MUOD  < ? > 
r 

3 6Q  CONTINUE 
r 

1 F ( KCQUNT  ♦ GT-^-24  CO  TO  390 

ntniC7  = n, 

QLDCHZ(l)  = 3.  - - - - 

OLnCHZ(2)  = C. 

QL30D7  < 1 ) = 0. 

OL  noo  7(2)  - C. 

.■a,a....lJS.Z.,.l,=  .!*.NDJ3 

OLDOPZm  = c. 

362  CCNTINUF  

79r  SING  = CINC  ♦ .5*  ( XINC  ♦ 0L0IC7) 

TCOH(l)  = TCOHtl)  ♦ ... 5*  CXCDHi  1)  _*  QLDCHZUH 
T C OH  ( 2 ) = TCOH (2)  ♦ .5* ( XCO  H<  2 ) ♦ OLOCHZ(’)) 

T.MiiQ2..ni  SMunoin  * .sMxmxju  ^-otoaozn  > ) 

SM'IOOd)  = SHUOO  ( 2 ) ♦ . 5 * ( X 00  ( ? ) ♦ 0L00D7  ( 2 ) ) 

DO  392  I = 1 •»  NQ  R 

HOPINGd)  = OOPINCd;  ♦ .5M00PXCI)  ♦ OS.DOPZCI)) 
OLDOPZ(I)  = OOPXC I) 

CONTI NUr 

QLDICZ  - XINC 

OL^CHZd)  = XCOH(l) 

0L0CH7(2J  - XCOHF24  — - - 

01*007(1  ) = X00(1 ) 

OL  non  7 <2  i - X OD ( 2 I 


A -7 


| 


IF  U7  .f.E.NZHl  SO  1 0-  4-00 

Ip  ( 70  .LT  , TESTL  ) GO  TO  t*30 

ton  WRIIf  4«K^FI.-  — - 

410  FOnMAT(*l*t ’RANGE  CELL  NUMBER  * , 12  »5X, *C£LL  LE NGT  H*» 1PG9.2  , 

1 ?X,*JFT  PAOIIV/T 

00  404  1=1, NOB 

POE»lCL»  = 0Qfl-IMC(I)TABiAT4.*gI 

ORSOOPfl)  = 0 BFT N ( OOP  INC  1 1 ) ) 

TOriPUT  = TDOPII)  ♦ OOcINC<  I> 

FI  = T 

VELOCP  = PI^OOPVEL 
coooOP  = VEL00P*0ME  G/PI /C 

M8W4MW4  I Wf  LOOP  ,ffiSQC-P.,U4PlNC<  I)  ,OaSOOP<I> 

4C2  FooMAT(*"OOPPLrR  BIN  * , T 2 , 5X , *OOPPLER  *, 

1 'VFLOT.ITV  <M/S€C»  = *, 1P&1 Qw3,5X,*FR€QUENGV  SHIFT*, 

g— 1 *4HZ+ » * , G 1 0 . ?/T  1 E i *~S HHr*y-t X , *GH*  < T10  , l-P-110^3/ 

7 IX  ,*OBSM*,T10,1PE1 0. 3/) 

404  fONTINlF 
r 

SINC  * 5 1 NC  *4  « *PI * A PE  A 

OPS!  = OBFTN(SINC) 

900H«»4 , *0 1*  ARCA**g*  i TCOHf  1 » *-»  E»  T€OH< »»•*-*»> 

OBSC=n9FTN(SCOH) 

S<F>  ^ AOQ*ABEAOQ*<SMUO0<  l>-*SHU©t>«i»  ♦ SMOOO  ( >t  *SHtWo1  i M /*! 
OB rOO  = OFF TN ( SOD  I 
S7^*STHC  ♦ SCOH  * c00 
OBSMrOPFTNtSIG) 


r SUM  C^OSS  SECTIONS  F° OB  EACH  range  cell 

r 

TSFQH  = ISCOH  ♦ SCCH 
T cOO  * T SOD  ♦ 500 
T3TNC  = ♦ SYNC 

TBTf,  - T 6K,  » 510 

r poj*t  t/)TAL  CROSS  SECTION  FO®  S‘THTGE  "Cftt 

r. 

HP ITCIR, 415 > K«K 

4 1 c FOoMflT(*QTOTAL  CROSS  SECTIONS  F OR  RANGE  CELL  *,I2) 

HPITF.  <4,430  SlNG»SCOH,SOOfSlFjTOOSIfOOSOT JOSOPjfrOSH 

4’0  F0°MATtT15,*SINC*,T25,*SC0H*,T?5,*S00* ,r4S,*RCS*/ 
*iX,»SH*TT*5,lP4€l«,  3/-H<T*OOSH*f  m»l<N,£l<»v3A>  — 


r 


KK-K  -■  MMX.  * t — — 

TF^TL  - TFSTL  ♦ FL 

SfrMO--  0. 

TCOHd,  = 0. 

TCOH<?T  -»  <?  . - - — 

SMUODf!)  = 0. 

-SHWWW  -*-<=.-  - 

no  <*25  JJ  = 1 tN0° 

OOP-T-NG-FtM-T--*-  -9-. 

u?c  CONTINUE  • 

KOOWF  « « - 

r 

r;ONTlW!c 

r 

r r.n!NT  TOTAL  C,(* 050  OCOTIONO  CELL'S 

r 

wpitp<f,  *??> 

4’?  FO^MAMM  TOTAL  OOPPlfP  CPOCS  FACTIONS  FOP  ALL  RANGE  CELLS*/) 
00  >-*5  T*i,NO» 
rtQr00Ptl ) = 0«rTN  ( T OOP  ( I)  ) 

F4— { 

y/rL0Cc  = f:»oopvfl 

FROOGP  « V'LOGT^GG-G /P»y?^-9<i€^»fl 

WRTTf  (:,,*?»>  I,VELOOP,FPOOOP,TrOP(I)  •GnTOOJ(I) 

rooMAT<*9IV}PPtE^  GIN  *,  ! ?f  ^ * dOQF^.f  R V=tOO!f  Y tM/SFCt  - *, 

1 1 n 0 1 0 . 3 t t'fi*FRFOUFMCY  SHIFT  (HZ)  r*,G10.J/ 

?d^f*-FTW0»/!KT»5H*Tm,lPfl&t  5/1 

CONTINUE 

*■  GofT^TOCOH) 

OOIGO  » OPCTW4TSOG) 

OPTING  = 03rfN(TGI  NC) 

00  7 err.  « OOFTHF  TO  JG) 

W°TTr  ,uu C V 

4vr  co=MAT(*0 TOTAL  C°0S5  S-GTJGNS  F0»  AtL  »*NG‘~Grtt5*f 

MOTTr  f.  ti.er  ) Tei^c,  TSCOH,  T<-0D.TSIG,03T  INCtOBTCOH,  D9T00,0?TSIG 

i—Z£. — co^MA  T < 1 1 R T * T S (*4G*  t T^Gt * T SF-O***  t T T *T  SGG*f  T4yj * TRGG*/ 

liY,*SH*,TlC,lPHElC.:i/lY,*0P$‘1*,Tiw,lD1,£ll.  1/1 

<" 

r.P  TO  ICO 

FN.n  — • 


A-9 


e c 


i 


I 


i 

I 


i 


r * 


C.  . 


r. 


SU^ROUT INF  RAY  (XC, ZC »OMEG» R JET tSTEP» IPRINT , OOPVEL »NOB,L AM30A> 
•,PLF«1D0,25), IATT.IOEX,IOVD,PLVtl03,25l ,NNR(1C0),  IPUFF, 

2 ,ppif<  ioor2Si-r^g-tAU4c;i^&)  ,ga.R u. 

COMMON/RAY/STILOA,TTIinA(2>  ,FMU00(2)  ,DOPSTR(5D)  ,PLR0VD(1JC  ) , GAr.M  A 
COW*ON/ANCLf/TH£TATTUEIADrT-Hi-fAOVTM€TG4^COGRHfTSTNTH£TTf++ETA 
DIONS’ ON  AT(3),ATTN(3), PH( 3), PHSF (3), PSCi:C) 


-RE  AC  J^WSTW*ERTTl--rNfWTT2,rNEWATTrN*WPH,NOOTK,NOOT¥tNOPMEX,NSTL  , 

1 ISManfl 

aAXft-CA3^8W*0^X  , FNC/3-lA.E*>-U/r  0-1-/*  . 14lS»2fc4>/ 

0 A" A CUT0Fr/-2G./ 


r DETERMINE  "-MIN  AND  SMAX  INTEGRATION  LIMITS 
r.  PLUME  DIVIOEO  INTO  N7M1  DISCS, 

L CC°  TN  I INCREMENTS  DISCS. 

.C-.S.AY.  I-E ?. I E-0 - f .Ofi- im.2 ftS EC T I C L UlfH.QISC 


77  C = 7G  - 

TPtTHPTA.GT.THFTAOl  GC  TO  152 
17"  - 

IF(THRTfln.GE.l.»  ZZO  = ZC/TTHETA 

_ l52-.blf.ACE  ---3 

SMTn=i.E6 
SMAXi-i,£6 
NZM1  = N7  - 1 
CO  lie  I - 1 » NZM 1 
r 

Q DISC  RAOXUS  - MAXIMUM  PLUME  RACIUS  AT.  AXIAL  STATION  I OR I4j. 

r 

R=AMAXi  tPLR  (I ♦ NRI  II  ) ,PLR<  I*  1,  NRU*i  T)  ) 

PS  f T I :: 

r.  . 

C TrrT,  -'fly  MTSSCE  DISC  IF  X?  G Pr  ATf  K THAN  P 

C 

Ic  (XC.GT.O)  SO  TO  11C 

IF  ( AASISX  NT MCI  ,LT«,vl>  -40  TO  i« 


01  STANCE  FROM  PAY  INTERSECTION  NITM  X»7-PC4N€  AT  <Xd,ZZO> 

TNTr PSR^TIONS  WITH  INFINITELY  LONG  CYL INOAR  WITH  RADIUS  EQUAL  TO 


r otsc  rao ius  a 


S=SQRT<R**2-X0**2»/SINTR£ 
00  ICC  J=  1 , 2 
IF  (J.E0.21  S=-S 


A-10 


i 


r 


r.  J^LOOP-  -OFTreMIMES  Z ^ALUFS  OF  POTH  INTERSECTIONS  OP  RAT  WITH 

r TNFINITELY  LONG  CYLINOAR.  A TEG  T IS  MAGE  TO  DETERMINE  WHETHER  OP 
r-  NOT  TH-r  INTERSECTION  Z-VA-UJt  IS  WITHIN  THE  -OISO  Z VAT.Or  H.IWITS 
r 

Z=77?»S*COSTHF 

IF  <(PLZ(I+1>-Z)M7-PL7(I> ) .LT.C.)  GO  TO  1)0 

SMAX=  AMAX  1 ( SNA  X,  S > 

SMTN=AMIN1(SMTN,S ) 

10?  CONTINUE 

144  ■■O-QWT-y-NHC 

r 

r t w-  RAY  MAY  INTERSECT  A DISC  X,  V PLANE  FfrCE  frf-SHflr  OR  SMAX. 
r Tur  j L00P  Or  TfDMINES  p AY  INTERSECTIONS  WITH  THE  X,Y  PLANt  AT  EACH 
G AXTAL  GTATton,  THE  INTERS POTION  POINT  (AT  »AOItt)S  R1  IS  TESTED  TO 
- WMF  r h ” - 0»  NOT  IT  LIES  WITHIN  the  ^ISC  (OF  RADIUS  RS(I)> 


IE ( APS(COSTHE) .LT . . Gl)  GO  to  150 
RF(HZ»-RS<NZM1) 
n0  lU1'  T=1,NZ 
S*(*>L7(IT-77eH  /COSTHF 
°=GQOT(X0*X?+S*S*SINTH- *SINTHE) 

— {4-  ho.t-TtCE  *B  » CO  TO  1-C- 

ihi.l  .1)  nface  = l 

SHAX*Amax1<SHAX,S) 

SMIN= AMIN! (SMIN.S) 

!<•*  CONTINUE 
lEf’  C 0 NT  I NUr 

— 4 

r INTEGRATION  LC^P 
r 

novr  - icvr 

IF (THETA ,LF .THFTA6)  IIOVO  * } 
ron=  FNC*CMfG*OMEG*SINTHE*rINTHE 
OH  A SE  C » ■ -g . *«HE G *Pi>€  T « < E H 1*4  ♦ z 74*005 THE 4 / C 

I OUT  = C 

PEL A*  - ° IXLAHR9A 

TR  r ", 

T?  « 0. 

9OLO  = SOPT (X?*X0»70*7‘ > 

— OOSAfcO  — W 

I F ( ROLO  • GE..CT1)  COSAL  0 = 7. O/POLQ 

to<;  . -l , G*  SI  NTH?  *C  OS  ALP 

touts  - GTNTHE*SQRT (1 , -COSALP*COSALP) 

T 7 ^ - cnSTUF  

PHASES  = PHASEC 

OQ  Hi  I-l.NOP 

DOHSTPII)  = 0. 


A-ll 


o c 


I 


I 

I 


I 


I 


* 

i 

# 


I 

I 


OS  = STfP'PJFT 

XK  « -U — 

J = 0 

HC-— -9 

ion  = - 

f — «__Q 

IF(SMIN.GT.SMAX)  RETURN 


If?  I r 14-1 

« 

r PAY  COORDINATES 


C.IM1  a 4-4 

S - SHIN  ♦ FIM1*STEP 

Y = -1.»S*0INTHE 

R = SORT < X*X*Y*Y) 

Ig  a-MgTA-.4I  .THC.U04  M T-O- 1 (*, 

R = POLO 

IFTKLE.4T  00  TO  Ibb 
OS  = STrp*pjfT*TZ 

7 - ZOLD  *JZ*ST£0-  

R = POLO  *■  TP^STEP 

IflR.fcT^aq  £0.1-0-444 

R = A9"(P) 

TRS  = A8S(TPS> 

R = AP^P) 
l&F  ZOLO  = 7 
POLO  = R 

_C 

r INTER  INTERPOLATES  LAPP  INPUT  TO  DETERMINE  THE  ELECTRON  OENSJTY 

A NO  GOLLIS^ION  FREQUENCY  -A-T  T«f  POINT  -U,¥,7L  — 

CALL  INTER  < R*  2*€S*F:Sr¥£L~rLAN40A~rORAOF-Z  * GRAOF 7 » &RAOfcR^ONAOF N»IOOT  ♦ 

r 


r 


r INDEX  ''ETPPMINES  THE  INDEX  OF  REFRACTION  N s XlNO 

r 


CALL  !NOEX<ES,FS,CMEG»XINO,YIND»inSX) 


I YINO 
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-C-  GRAQIE-fiJ-  OF  IWW  -Of  - Pg-F-RACIXXlN  - 


RRAKET— =-  CMf**OM€0  * FS*FS 

JPf,  = i. -ES/FNC/BRAKET 

xaAaujJU  w = 

FN  = SOPT(APG) 

GRAONF  = G.  E/FNG7  BRAKfTyfN 

G9A0FF  = 2. 0*ES*FS/RRAKET 

GPAOXR  = GRAOWFMGE4GFF*GRAGFB-G*AO€R»  - 

G°A  0X7  = GCA0NFMCPA0FF  »GRA0FZ--GRAQEZ) 

FI--— f5*.f-S/0>»f&/0M€G 

F 2 = C.-ARG)M1.-A0G> 

F3  = SORT  I APG  * ARG  ♦ F2»F1> 

G°A0NP  = F2*F1*GRA0FP/FS 

GRAONR  = GRAONO  +•  GRAOXR*<ARG» < i.+Fl)  - FI) 
GRA0NO  = .25# (GRA0NR/F3  * GRAPXRJ/XINC 

GBAOM7  - S 

G R A 0 N 7 = G°A0NZ  ♦ G RACX 7* IA RG* ( 1 . *F1 ) - FI) 
GPAQN7  = »?F* (GRAONZ/F3  * GRABXZ)VXIMO 
GOftOW*  » GQRT  <Gft*PMR«6RAPN»»OftAONZ«GRftOMZ) 

r 

r UNT?  NOP«A*.  VECTOR  0 

r* 

TP  ■=  T»G 
T 7 = T?E: 

TPm I » TP-HT-S 

UFACTo  = TR*GR AONP  ♦ TZ*GPA0NZ 

l>o  * G»A«NR  - MF ACTR*TR 

UP  = UR/GRAPNM 

OZ  * GPAON7  - Of AC TP* T 7 

U 7 = U7/GPA0NM 


I 


r 


P UNIT  TANGENT  VE€T OR  T 

C 

FOTRS  = GPAONR'OR  ■+  ORAGNZ*OZ  — — 

cOTQS  = FOTDS/FN 

OPngR  r cotos»UR 

DTnS7  - POTOS’UZ 

OTOSP  * POT05*uphj — 

TPS  = TR  ♦ OTOSR*STEP 

TZF  *-  IZ  ♦ QTOSZ*SF£P 

TPHIS  = TPHT  ♦ OTOSP*STE° 

XHAfr-a  SORT  »TRS«TRS  ♦ TZS*T-ZS •»- S*4  PH  I 5 ) 

TRS  = tps/THAG 

-—  — TZ^  ^=-47S/T-HAG 

TPHTS  = TPHIS/TMAG 

- SVEN  I 


1 F8  IF (HOD(I « 2) .NE.O)  GO  TO  171? 
SIMP  * N, 

GO  TO  180 


r ooo  r 

£ 

l?p  SIMP  = 1. 


0 PPEVTOUS  I NT c GRAND  AOOEO  EVERY  OTHER  POINT  FOR  SIMPSON'S  DULE 
18"  TF(I.EO.KK)  GO  TO  19C 

Oinsn-  « G. 

OLOTT1  = 0. 

Ot.«TT?  ^ F.  — - 

QLnATT  = C. 

04.  riPW  - 5 . - - — 

GO  TO  ?01 

IRC*  KK  = K*  *■  ? 

r 

r.  SIMPSON'S  ®Ot-P — 

r 

ZO*  NFWATT  - ATT*OS*5INP/8, 

ATTrNS  = ATTENS  ♦ NEHATT  ♦ 010ATT 

04. PAT  T » NCWATT 

NrW°H  - Q*OS*SIMP/7. 

PHASES  •»  PHASES  + NgWPH  ♦ OIPPH  — 
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QLHPH  -.NEHPH 

A T T EN  - 4 TT  FNS 

PH  ASF  = -P-HASF5 

r 

_C-  SIMPSON  2 PULF - 

IFU.LI.4)  GO  TO  213 

I F ( MOC (1,2) .NE.O)  GO  TO  21f 
AT^PIS  = UTU)  + 3,*AT<24  * 

A T T CN  = ATAPTS  A T TN  (5 ) 

— - PHAP-TS--  4PW43I  * 3.*PW<2)  *■  3\»PH41>  «•  Q)*.3?5*0S 
omasE  = PH4PTS  ♦ PHSE  < 7 ) 


0 Ti-r  FI^ST  ANO  SECOND  POINTS  ARE  SPECIAL 
0 

2lr  TF(I.F').l)  ATTEN  = 0. 

IF(I.EO.l)  PHASE  = PHA5EC 

IFCI.F".'M  ATTEN  = OSMATC)  + ATT)/2. 

IF  tI->£0.  21. -CHASE  =■  .CS * XPH.U4— fc-OXy-Q. ■*  PHASES 

r 

r.  STOPS  LAST  THSgP  POINTS  — 

r 

IF(I.FO.l)  GO  TO  233 
I F ( I , FO . ? ) GO  TO  2Zr 

A T ( 3.) j,.  .A  l l.  21 

AT^F(T)  = 0 TTN ( 2) 

PHU)  - PH(?> 
pp^EC?)  = pHSE(2) 

22C  A T ( 2)  = A T ( 1) 

AT  ▼ M ( ? T r ATTN(l) 

pHf?l  - PH(i> 

PH  '*£'('>)  = PHSE(l) 

23C  AT ( 1)  r ATT 

ATtn<1>  - ATJEN 
PHtlJ  : 2 
PH<^F(1)  = PHASE 


r A TT  pNU  AT  T ON 

cryp  - 2 , ® A TT  EN 
EX  = l./EXP  (FEXP) 

Ipf TATT.EC.C)  FX=1. 

£ 

C TPPMINATICN  CONDITIONS  i OVtRDf NSF  SURFACE  CROSS  SECTION 
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Ir(IOO.EQ.l>  GO  TO  231 
IF  ill  OVD*  EQ  .C  ) GO  TO  2<*£ 

RAT  = rS/Cfc'rG 

ZNGB-jl-FOO*  (1 . + KtT»BnUtU  »gA34 

IF(ES.LT. 7N0D)  GO  TO  2-2 

inn  = i 

231  I F f N'r  ACE . EO . 1 ) GO  TO  250 
no  232  N^2«N7 
I«r(Z.Lr.PL7(N)  ) GO  TO  234 

22?  CONTI MU£ 

2T<*  IF ( N . G T , N 7 ) N = N7 

TANPSI  = FLR0Vr<N4  - PIROVPIN-l) 

TANP-S-f  » TANPSI/ (PL  7 Pi?-(N-1-H 

COSPSI  = l./SQRT( 1. ♦! ANPSI*TANPSI> 

SIMPRI  = SQ°T (l.-COSPSI*COSPSI) 

XI  = 1. 

TF< IRiFF. PO.O)  GO  TO  23S 

PfT  = TANFSI*(Z-PL7(N)»  *■  FL~OVO(N> 

U - ,fc*pjr. t 

IP(PCO.LE.H)  GC  TO  235 

XII  =•  .A*OHEG*»J£T*SINTP?/r 
XIT  = XII*XII 

~W  = M . ♦ »10**I  I /POD/ROD)  *£XP (— 1 . *Xl  I > 

XI  = cORT(YI) 

Y-STEP/2, 

RTMqifl  = Y/SQRT(X*X*-Y*Y) 

NOOTK  r -1  .*SINTHF*COSPSI*S  It^BTft  - COS  THf*5INPSI 
NOOTK  = A RR (NCOTK ) 

NOPTY  - COSPSI*SI N°TA 
P A T T 0 = NOOTK/NOOTY 

— Of  NO*  » ( 1 1 ♦*  INO+*<t"> +«  I*OV  +-  *IN0«V!NP 

T'-YMAG  = 1.  - XI  NO  * XT  NO  - YlNO*YIND 
NO°Mfy  - CUFG*°AT IC*PX*XT/T1 

FMUOOC.I  = NOR^EX*(OEXMAG*COS(PHAR:)  ?.*YINO*SIN(PHASE > > /0ENOM 
P“UOP<">  ■=  NORH£X*<2.*YINP*COS(PMASE)  - Of  XMAG*SI  N(PHASS ) t XCfNOfl 

GO  TO  "50 


r LA'-t  i MUr-r  O'-  ODD  FOP  SIMPSON'S  RULE  INTEGRATION 

242  IF  CIOUT.  pp.  1)  GO  TO  25.' 

I-MThFTA.GT  « THE  T A 01  GC  TO  ?5=r 
OP'X  = 12 ,’ALOGIC (AMAX1 (EX, l.E-1 ') ) 

IPMWCY^-f-t-CUTOfr  » PO  TO  050 

OTIS*  = Q. 

IF<T HE TA,tE»T HE T A 0 ) OTOSH  » SORT  ( CTOSR*OTOSP«-OT05  Z*GT05Z'*- 
l o 7nsP*oT0EP) 

GO  TO  EGG  ------  — 

25n  IP (MOO(It?) .NE.0)  NS  = I 

2-Cg  OONTINOF 
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w<=TTr  (rit?qc  ) S.£X»5TinA,TTlL0fl(ll,  TTILOA(?»  , 

i PMuonmtrMuoo<e»,ooFST=<wc8) 

->q|'  r0"M£  T (1H  , 1°9G15 . 7 ) 
r 

’0'  go  to  ibc 

r* 

r "NO  OF  INT"GRATION  LOOP 
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IF  (I«GT.l)  GO  TO  310 

STILOA  - C._  - - - 

TTILOA(l)  = C. 

III4.OA.UA--.  4U 

31?  CONTINUE 

R.rjUPN  — - - - 

£NC 

SUnROUTINE  PR0F( I,P,F,F,VP,LAMODA,GRA0ER,GRADFR,I0UT) 

r ocOF  IF  CAILEO  OY  INTER  TO  PERFORM  iNT  ERPOL  A TI ON  IN  R. 

C THE  INTERPOLATION  IS  LINEAR  --FO®  INF  GOLLISSION  F^eOUENOY  A NO 
r AXIAL  VELOCITY,  BUT  QUAORATIC  IN  TH E ELECTP ON  TENSITY 

c_- 

COMMON/PLUME/NZ,PL? (ICC) , NR ( 1C  C ) ,PL R ( 1 C 3 , 25)  , PL E( 1 ' 0 , 2 5) 
*,PLF<  100r2S>,  I ATT  , I OE  X*  IOVO  *PL  V < LOG  ,25)  T*N-mM ) , IPLFA, 

2 P°LE  <10° ,25) ,PPLF< ICC ,25) ,P=LV(1C~ ,25) ,3PlR(i: 0, 25) 

REAL  LAMOOA 

IP  (R.G’.PL^ 'T, 1) ) GO  TO  1C 

E - GLF  < I , 1) 

F = PLr ( I « 1 ) 

VR  = PL V ( I, 1 ) 

GPAOER  = C. 

GRA3FR  = C, 
pc  TURN 

— 13  JMAX  - NR  (I) 

r 

r.  TNI S LOOP  OE T FRMI NFS  PL R I I*J)  SUCN  THAT  

r PLR(I, J-l) .LT.R.LT.PLP(I,J) 

r 

CO  ICG  J = ? , JM A X 

LF  (B.LE  . PLB ( I ,J)  ) PC  IQ  113 

10C  GONTINIJP 

IF (R,L~ • PLR  < I , JNA  X) * , 1 ) IOUT  =•  1 
CELP?  = P - PLR(I.JMAX) 

A - 1,/LAMflCA 
AO  = A*P"L0? 

EM  - PL£  t I,  JrtAXJ 

°A  T IOF  r 1, 

RATIOV  = 1* 

IFfEM.CE.l.)  RAT  I OF  r PLF C I , JMA X ) /EM 
I F (EM,Gf * 1. ) RAT I OV  = RLVA I , JNAXT/EN 
F - FM*E XP( -1  * • AO ) 

IFtE.lT.  1 .)£■!. 

F r ° A TI OF* E 

VR  - RATI CV*E — — — 

GR  AOcP  = -1  ,*E*A 

GRAOFR  = GR  AOc-R  *R  A T-1QF — 

RF  TUPN 
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^ l i j Jiwn  ’ ’vrvunTvn — inro  hairl — rrruwTT# 

r THF  INTERPOLATION  IN  R IS  DONE  9Y  CALLING  PROF, 
r intcr  RCRfOR**;  LINEAR  INTERPOL  A TIDH  IN-?. 
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*,PLF<  100 ,25),  IATT,  IOEX,  IOVO ,PL V < 10 3 , 25 > ,NNR  ( IOC  ) , IRUFF, 


REAL  LAM90A 

00  100  1-1,07  — — 

IF(PlZ<T>-7)  100,110,110 

104  00HTI-N0F 

TOUT  = 1 


II  = N7-1 

0AL-L  ^^OF^I^rR-rEvFrVT-LAHOOAiOPAOEPtGRAO^fl^tf^ 
RATIOF  =1. 

RAII0V  -=  I, 

IF(E.Gc.l .)  RATIOF  = F/E 


IF ( 12 ,LE • 1 l GO  TO  115 
0AU.  FCOFTI1,0t€1  y^i, 

DEL72  = Z - PL7 ( 1 2) 

-A—*  -1,/H*R0A 

AO  = A *0 EL  7 2 


IF(E.LT.l.)  E = 1. 

F * RAT IOF*f 
\l  = PAT  10  (/*  E 

0OA00?  *L,«E  »A 

GRAOFZ  = GRA0E2*RATI0F 


11"  I ?= I 

CALL  P°OF (I2,R,E«F, V, LAMBDA ,GR AOER, GRADER, I OUT ) 
RA  T I0F  * 1, 

RAT IOV  “ x. 


IFte.G'.l.)  PATIOV  = V/E 
-IF  T I g-GT,  11  GO  TO  120 
115  OECAY  = E YP (Z/L AMB0A1 

E = E»OE0AV  

F = RATIOF*E 


GPAOEZ  = E/LAMPOA 

G®  AOc  I RATIOF^GRAOEZ- - 

Rp  TURN 

120  CALL  < 1 1,  R , EItFI , VR1 

OELZ  = PL Z < 12 1 - PLZ(Il) 


E = ( El  *OEL  1 ♦ EMi.-OElIM 
F - TF«-*-Oftl  ♦ F*<1,-B£L1M 


READ  (1,1CC)  TITLE 


RFAP  (1,110)  BUFF 

IF  (EOF ( 1 ) . NE • 0 « ) GO  TO  190  

IF  (ARS(0UFF<8)> .GF.1.E30)  GO  TO  160 

IF  U-l>  170,120,130 

l?"  IF  (PUFF(l)  .NE.O. ) GO  TO  190 
r_n  rn  i un 


130  IF((«UFF(1) .LE.OIO) .OR. (BUFF(l) .GT.10.) ) GO  TO  193 

1 GO  OLO=RUFF  ( 1)  

PPlR(I,J)=OLD 

IF (RUFF (2J.LJ.1.Q)  BUFF  124  —=  1.0  

Ic( BUFF(3).LT ,1*0)  BUFF ( 3 ) = 1.0 
TFJJIR5L1B11FB1AL1  flilF  F i A 1 - 1 . n 


PPLV(I,J)  = BUFF ( 8) *. 3048 

PPLE(I.J»=BUFF(2)  __ 

150  PPLF( I,J)=9UFF(3) 

J-27  __  

GO  TO  193 

if>a  im&ouui 

IF  (BUFF (8) .L T . 0 . ) GO  *0  180 

170  CONTINUE  - - 

1 = 5 

GO  Ta  190  


I 
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RETURN 

-19C-  IE5=1 

WR ITF  (f- , 200  ) I,J 
20f.  FORMAT  </iH4^*P4PUI 
RETURN 

POE) 


SUBROUTINE  INUEX(  ZNE,  ZNU,-OMfS,X,Y,rOEV) 
If(IOfX»  20, 1C, 20 
10  X - 1 , 

Y = C. 



?C  RAT  = ZNU/OMEG 

YY-ZNf/0^er,/3.  l<*E-13/eM£G 
YY  = YY/(1 . *PAT**2) 

XX-i , -YY 
YY=RAT*YY 

X*  fSQ»T-«-XX»»3*  YV»»3>  > XX  ) /g. 

X=SG°T(X| 

Y « ,5MSOPT<XX*XX+YV*Y¥»  - XX> 

Y = SOPT(A«S(Y»» 

RFTUOM 

rNn 


SUaROUT  rH£  SUNOTZNF^ZWtG AMMA T G3E T A tOM£G,C4- AM,&SiG,AIT,  — 
*OTAU,ETA, XIN0,YIN01 

DATA  Pl/4»tMS9g*&W - 

ZK-0MCG/,.998E8 

GAMSQ  ■= — GJUiMAAEAiUlA 

CK0LMrl5.6*CLAM**3/<i.^|2.* XIN0*ZK*CLAM)*>2) *• 1.8 33 

R£7Wf x7.817A€-9*ZNF  

OTAU=pcZNt/SQRT(l,«-  (ZNU/OMrG»  **2> 

OSI&=OTAU**Z«F-KOU«*GAKSQ  — 

FT A=?, *XIND*Z< 

ATI-ZKYYIAO»PltQSIGMl.»Ca£TA) 

R*-  tijon 
2 NO. 
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I 


I 


I 


****** ******  » **********  ******  T. 


s<>oroi>ti*e  cyu>€R<iep,titlsi) 

C^MM0N/PLUME/N7,PLZ ( 1 uC ) * NR ( 1 C : ) * PI  R C 1 C 0 , 25)  , pLE  ( 1C  0 , 2 5*  , 

1 PLc(iCe,25),IATT,Il}£X,rOVO-,f>tV<l-ee,25),NNR<lfl-et,  IRDFF, 

2 PfLf-  (IOC  .25)  ,P?LF(  10 G ,25) , PpL V < ICO , 25 ) ,PPI_p<1G  3, 25) 

OI*"NST(W- ^ITLrl(  20  ) 


p PONSTANT  CYLINDER  TEST  CASE* 

C LrNGTH  = N7-1  JET  PAOIT 
p PAOIU5  = . 1*  <NRR-i > JET  RAOIT 
p ELErT3ON  DENSITY  = D (CH-E) 

-G COLL  I SSI  ON  FREQUENCY-  — F f-H-?-) 

p AXIAL  GAS  </r  LOCI  T Y = V (H/S^C) 

r 

TATA  N7  « NRR/3  *11/ 

OATA  A,F,V/l,*>«»,S.i*ll»l,/ 

r 

Pf-AtMi-r  !<?■»■  TITLED — 

IP  POPNAT(20  A«) 

PL  * N 7-1 

PP  = .1*CN»C-1) 

H»IT?<~,1P)  FL,FR,0,F,V 
1=  FOcHAT(*l CONSTANT  CYLINDER  PROPERTIES*/ 

i*  LENGTH  -=»  rl-Pf,lflT.T  .gX,«JCT-4AAH-»/ 

2*  PAOTUS  =*,1DG10.3,2X,*JFT  RADII*/ 

?•  "LPCTRON  D-N'ilTY  •*fiPGl-1,3, 2X,*CH-T*/ 

POLL  I SS  ION  EPEQUENCX  =*  , 1 PC,  1 C . 3 , 2 X ,*  H7*/ 

•**  AXIAL  GAS  VELOCITY  »* , IPG  10  . 3 , 2 X , *M/  SEC*// ) 
Ir  D = ' 


CP!  7C  !=1.N7 

F T HI  * 1-1 
3L  7 ( I ) = r!Hl 
NNC(I)  = N»C 
NPI  = NNR(T) 


OP  2C  J=1,NPI 
PJ**1  - J - ! 

PPL°(I.J)  - P JM 1 * , 1 
0PLV<!,J>  * Y 
poL  E ( T . J)  = D 

PflLS44W» -*  F 

2^  rONT  I N!*P 
CONI  I NUf 
pr  T(jp»j 
ikjn 
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, APPENDIX  B 

SAMPLE  OUTPUT 
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APPENDIX  C 


INPUT  DATA  SHEETS  FOR  TYPICAL  RUNS 


t 


APPENDIX  C 


INPUT  DATA  FOR  TYPICAL  RUNS 


SYSTEM  A 37  KFT 
CARD  A:  15  DEGREES  (20A4) 

CARD  B:  15.,  1.  32E9,  .4,  5.  E4,  15.,  .65,  300.,  0.,  15.,  1. 

(free  format) 

6 (deg),  f (Hz),  JET  RADIUS  (m),  RANGE  (m), 
RESOLUTION  (m),  DIVERGENCE  (deg),  DOPPLER 
VELOCITY  (m/sec),  DOPPLER  FREQUENCY  (Hz), 
LOW  ASPECT  ANGLE  (deg),  GAMMA 

CARDC:  .1,  .1,  .5  (free  format) 

FRACTS,  FRACTX,  FRACTZ 

CARD  D:  0,  0,  1,  1,  1,  0 (free  format) 

ITEST,  IPRINT,  IATT,  IDEX,  IOVD,  IRUFF 
CPU  TIME  ON  CDC  6600:  1850  sec 


C-l 


SYSTEM  A 5 KFT 


CARD  A: 
CARD  B: 
CARD  C: 
CARD  D: 


CARD  A: 
CARD  B: 
CARD  C: 
CARD  D: 


30  DEGREES 

30.,  1.  32E+09,  .4,  1.  E+04,  5.,  . 65,  150.,  0.,  15.,  1. 

. 1,  . 1,  . 5 

0,  0,  1,  1,  1,  0 

CPU  TIME  ON  CDC  6600:  235  sec 
SYSTEM  D 23.  9 KM 

138  DEGREES 

138.,  1.  32E+09,  .621792,  3.  05 E+04,  15.,  .9,  340.,  0.,  15.,  1 
. 1,  .1,  .5 

0,  0.  1,  1,  1,  0 

CPU  TIME  ON  CDC  6600:  108  sec 
CPU  TIME  ON  CDC  6400:  344.  8 sec 


C-2 


